279 



VI. On the Law of Distribution of Molecular Velocities, and on the Theory 
of Viscosity and Thermal Conduction, in a Non-uniform Simple 
Monatomic Gas. 

By S. Chapman, M.A., D.Sc, Fellow and Lecturer of Trinity College, Cambridge, 

Communicated by Sir Joseph Labmor, F.R.S. 

Eeceived October 5,— Eead November 18, — Revised December 1, 1915. 



Contents. 

Page. 

§ 1. Introduction 279 

§ 2. Definition and preliminary consideration of the problem 281 

§ 3. The equation of transfer of molecular properties 285 

§ 4. The effect of molecular encounters 289 

§ 5. The general expression for AQ X 295 

§6. The form of the function F (u, v, w) 299 

§7. The calculation of AQ L 301 

§8. The expressions for the coefficients in the velocity- distribution functions ...... 311 

§ 9. Consideration of particular molecular models 320 

§ 10. Numerical calculations for particular molecular models ...... 327 

§11. Viscosity and thermal conduction 334 

§ 12. Comparison of the theory with experimental data 341 

§ 1. Introduction. 

The kinetic theory of gases can be developed accurately only after the distribution of 
the molecular velocities has been determined. This was done by Maxwell # in the 
case of a uniform gas, and by means of his well-known law of distribution the pressure 
and temperature can be precisely expressed in terms of the molecular data. His law 
does not suffice, however, for the investigation of diffusion, viscosity, or thermal 
conduction, since these occur only when the gas is not uniform in composition, mean 
velocity, or energy. An accurate theory of these phenomena must be based on the 
evaluation of the modified velocity-distribution function, a task which for many 
decades has constituted one of the classical unsolved problems of the kinetic theory. 

* Maxwell, 'Scientific Papers/ L, p. 377, II., p. 23. The proofs were unsatisfactory, and have been 
improved by Boltzmann, Jeans, and others. 
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In one special case, as Maxwell found, the actual determination of this function 
proves to be unnecessary for the purpose mentioned ; this is the case of a gas composed 
of molecules which are point centres of force varying inversely as the fifth power of 
the distance. The reasons for the peculiarity in this instance are analytical a,nd not 
physical, and unfortunately for the simplicity of the mathematical theory of gases, 
Maxwell's results # for such a gas do not accord with the observed data of actual 
gases. This particular molecular model is therefore interesting chiefly on theoretical 
grounds, and it is important to develop the theory for molecules of other types, 
which may better represent the behaviour of real molecules. 

Until recently no progress had been made towards the determination of the velocity- 
distribution function for a non-uniform gas, beyond a theorem by BoltzmannJ" who 
proved that the function must satisfy a certain integral equation. In 1911, EnskogJ 
applied the method of solution by series to this equation ; he determined the form of 
the function, but without evaluating its coefficients, and his numerical approximations 
proved far from satisfactory. In 1912, Hilbert§ showed that if the molecules of the 
gas are rigid elastic spheres, Boltzmann's equation may be transformed into a linear 
orthogonal integral equation of the second kind with a symmetrical kernel, and 
deduced the existence of a unique solution. Lunn|| and PidduckIF have since removed 
Hilbert's restriction to a special type of molecule, and by means of the transformed 
equation Piddtjck has worked out a numerical solution of a special problem . on 
diffusion. These researches are of much importance and interest, especially from the 
logical standpoint of the pure mathematician. The use of Bolt'zmann's equation, 
howevever, does not appear to be the best method of actually determining the formal 
solution ; thus Piddtjck states that the symmetrical kernel of the transformed 
equation shows no special properties in the case of Maxwellian molecules, and in the 
numerical solution it appears to be necessary to repeat all the calculations, which are 
very laborious, in every special case which is worked out. 

In 1911, by the assumption of a simple form for the velocity-distribution function, 
I endeavoured to extend Maxwell's accurate theory of a gas to molecules of the most 
general kind compatible with spherical symmetry. ## Subsequent acquaintance with 
Enskog's work convinced me of the approximate nature of my results, and during the 
last few years I have given much thought to the determination of the general velocity- 
distribution function. By a method which is quite distinct from that based on 

* Maxwell, 'Scientific Papers/ II., p. 23. Molecules which are point centres of force varying 
inversely as the fifth power of the distance will, for the sake of brevity, be referred to as Maxwellian 
molecules. 

t Boltzmann, ' Yorlesungen iiber Gastheorie,' L, p. 114. 

| Enskog, ' Physikalische Zeitschrift/ XII., 58, 1911. 

§ Hilbert, ' Math. Annalen/ 1912, or « Linearen Integralgleichungen ' (Teubner), 1912. 

|| Lunn, 'Bull. Amer. Math. Soc.,' 19, p. 455, 1913. 

U Pidduok, 4 Proc. Lond. Math. Soc.,' (2), 15, p. 89, 1915 ; cf. p, 95 for the statement quoted, 
** Chapman, 'Phil. Trans./ A, vol 211, p, 433, 1911, 
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Boltzmann's equation, viz., by the use of the aggregate of the equations of transfer 
for certain infinite sequences of functions of the molecular velocities, an expression for 
the velocity-distribution function similar to that found by Enskog can be obtained, 
and general formulae for the coefficients can be determined. The present paper contains 
the solution for a gas in which the mean velocity and the temperature vary from 
point to point, the results being worked out at all completely only for the case of a 
simple gas ; in a later paper I hope to give the solution in the most general terms, 
so as to yield a complete theory of viscosity, thermal conduction, and diffusion in a 
composite gas formed of two kinds of spherically symmetrical molecules of any type. 

The formulae obtained by the present method lend themselves to numerical calcula- 
tion, and are found to converge rapidly. The results for any particular molecular 
model can be calculated to any desired degree of accuracy ; in this paper three special 
types of molecule have been considered, viz., point centres of force varying inversely 
as the n th power of the distance, rigid elastic spheres, and rigid elastic attracting 
spheres. It is found that, for such molecules, the errors in the approximate formulae 
for viscosity and thermal conduction which were given in my first paper do not exceed 
two or three per cent, at most. The detailed numerical results, and comparison with 
observed data, are riven in §§ 10-12. 

§ 2. Definition and Preliminary Consideration oe the Problem. 

The Nature of the Gas. 

§ 2 (A). The gas contemplated in our calculations is monatomic and nearly perfect, 
"monatomic" implying nothing more than spherical symmetry of the molecules, 
while " nearly perfect " denotes a certain state as regards density and temperature ; 
this state is such that the molecular paths are sensibly rectilinear for the majority 
of the time, being altered by mutual encounters, the duration of which is a very 
small fraction of the average interval between two encounters. In these circum- 
stances the number and effect of encounters in which more than two molecules are 
simultaneously engaged is negligible in comparison with the number and effect 
of binary encounters. 

The gas is supposed to be acted upon by external forces, and the variations of these 
forces, and of the density, mean velocity, and temperature of the gas, with regard 
to space and time, are small quantities of the first order at most. In the present 
paper the density of the gas is supposed such that the mean length of path of 
a molecule between collisions is small compared with the scale of the space-variation 
of the above quantities ; the modifications .of the theory in the case of highly rarefied 
gases, where the mean free path becomes large, will be dealt with in a future 
paper. As we are not interested in the mass motion or acceleration of the gas 
as a whole, but only in the small variations with regard to space and time, it is 
convenient to imagine that, by the addition of a suitable uniform motion and field 

2 q 2 
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of force to the whole gas, the mean velocity and acceleration at the particular point 
and time under consideration are reduced to zero, the velocity and acceleration at 
other points throughout the gas being small, though not actually zero. 

Notation, 

§ 2 (B) We shall denote the mass of a molecule by m, the number of molecules per 
unit volume at the point (x 9 y 9 z) by v 9 the components of external force acting on a 
molecule at (x 9 y 9 z) by (X, Y, Z), the components of the velocity of a typical molecule 
by (u 9 v 9 w) 9 and the components of the mean velocity of the gas at the point (x, y 9 z) 
by (u 09 v , w ). The vector difference between the velocity of a typical molecule and 
the mean velocity (u , v 09 io ) will be called the peculiar velocity of the molecule ; we 
shall denote its components by (U, V, W), so that 

(1) U = u — u 09 V = v—v 09 W = w— w . 

The Distribution of Velocities. 

§ 2 (C) The distribution of the molecular velocities may be specified by (u 09 v 09 w ) 
together with a function /(U, V, W), called the velocity-distribution function, 
which is defined by the following property : the number of molecules contained 
within a volume-element dx dy dz about the point (x 9 y 9 z) which possess peculiar 
velocities whose three components lie respectively between (U, V, W) and (U+cZU, 
V + dV, W + dW) is 

(2) v/(U, V, \N) dU dV d\N dx dy dz. 

Besides being a function of U, V, W, f will depend on the mass m 9 the absolute 
temperature T and its space derivatives at the point (x 9 y, z), and on the space 
derivatives of (u^ v Q , w Q ), but not on the absolute magnitudes of the latter : for we 
may evidently impart an arbitrary additional velocity (id, v\ w f ) } to the whole mass 
of gas without affecting the distribution of the peculiar velocities of the molecules 
at any point. It is therefore legitimate, and it will prove convenient, to suppose 
that, at the actual point under consideration, u () = v — iv Q = ; where u 09 v 09 w 
occur in any expression which has to be differentiated, however, they must not be 
made equal to zero till after the differentiation has been performed. 

In consequence of the definition of/ and of U, V, W, / must satisfy the following 
equations : — - 

(3) |/(U, V,'W)dUdVdW = l, 

(4) U/(U, V, \N)dVdVdW = V/(U, V, W) dU dV d\N 



w/ (u, v, w) <2u dy dw = 0. 
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If Q denotes any function of the velocity components (u, v, w) of a typical molecule, 
while Q denotes its mean vahie at the point (x, y, z), we have 



(5) Q = Q/(U, V, W)d\JdVd\N, 



in which, for purposes of integration, Q would be expressed in terms of u -\-\J, w + V, 
w +\N by (2). In the integrals (3) to (5), and elsewhere throughout the paper, 
integrations with respect to the velocity components are understood to be taken 
over all values of the variables, from — 00 to + °° . 

The equations (4) may, in the notation just introduced, be expressed as follows : — 

(6) = V = W = 0. 

The Velocity-distribution Function for a Uniform Gas. 

§ 2 (D) When the gas is uniform, all the derivatives of T and of (u , v , w ) are zero, 
and f must depend only on m, T, and (ll, V, W). It has, in fact, been shown by 
Maxwell and others* that 

f*j\ f ^ (hm\ 2 e _/ m (u 2 +v 2 +w 2 ) 

where 

< 8) k = RT > 

and R is the universal gas constant in the characteristic equation of a gas : 

(9) p = Ei/T. 

The Distribution Function for a Non-uniform Gas. 

§ 2 (E) When the gas is slightly non-uniform, f will differ slightly from the value 
given by (7), which we shall denote by f : we may therefore write 

(10) /(U, V, W) =/ .(U, V, W) {1 + F(U, V, W)} = ^JV Aw(u ' +v,+w ^ {1 + F(U, V, W)}. 

The function F will be of the same order of magnitude as the variations of 
temperature and velocity in the gas ; these space derivatives we shall regard as being 
of the first order, and as we shall neglect second order quantities throughout our 
work, no products of derivatives will occur in F. Hence, since F vanishes when the 
variations in the gas are zero, it must be a linear function of the space derivatives 
of T and (u , v ,-w ), with no term independent of these derivatives. The coefficients 
will be functions of m, T, and U, V, W. 

Clearly the form of F cannot depend upon any special choice of axes of reference 
(these are throughout taken to be mutually perpendicular), so that F is an invariant 

* Of. Jeans' * Dynamical Theory of Gases/ 
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with respect to any orthogonal transformation of the co-ordinate axes. This places 
some restriction upon the possible modes of occurrence in F of (ll, V, W) and of the 
space derivatives of T and (u Q> v 09 w ), though not, of course, on the scalar quantities 
m and T. It is easy to see that the most general invariant function of the quantities 
involved in F must be compounded of the following elementary invariants : — 

(11) C 2 EE U 2 + V 2 -fW 2 , 

(12) (13) 's = ^ + ^ + ^° v 2 T^f-^ + ~^ + ^T 
{ h { } - dx + dy dz > ~" W + dy 2 ^ dz 2 ) ' 

(14) DT = (u 4- + V #- + W %-) T, 

V ox oy dz) 

(15) S' = U 2 ^+V 2 ^ + W^ + Vw(^ + ^Vwu(^+^)+Uvf^ + ^) ( 

dx oy dz \dy dz I \oz dx J \dx oy J- 

together with derivatives of the last four expressions formed by operating on them 
any number of times by the invariant differential operators V 2 and D, in the notation 
of (13) and (14). 

[January 15, 1916. — Except in the case of highly rarefied gases, which were 
expressly excluded in § 2 (A), only the derivatives of the first order actually occur in 
F, to the present degree of accuracy. The reasons for this will perhaps be more 
clearly apparent after reading § 11, but the following considerations will elucidate the 
point. Whatever derivatives are contained in F must (§11) appear either in the 
equation of pressure or the equation of energy, so that, if the ordinary equations of 
viscosity and thermal conduction are to hold good, only the first-order space deriva- 
tives of temperature and mean velocity can be present ; otherwise the ordinary 
coefficients of viscosity and conduction do not exist. In actual gases at normal 
densities the ordinary equations are shown by experiment to be valid ; they fail, 
however, in highly rarefied gases because the terms in F which contain second-order 
differentials of T, u , v 09 w are in this case comparable with those containing derivatives 
of the first order, as will be seen in detail in the future paper mentioned in § 2 (A). 
The coefficients of the first and second order derivatives respectively contain (x/Z) and 
(x/Z) 2 , where X is the mean free path of a molecule and I is comparable with the scale 
of length within which the temperature and mean velocity vary appreciably ; except 
in rarefied gases (x/Z) 2 can be neglected in comparison with (x/Z). The same inferences 
can be made also (cf. § 6) from the equations of transfer of § 3. 

For the present paper it is therefore sufficient (and it is convenient) to write down 
the following form of F forthwith : — 

(16) P = (u|^ + v|^ + W^P 1 (CP)+SP a (C«)+S / P3(C 3 ), 
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this being the only way in which the first-order derivatives can occur in F, in order 
that F may be an invariant. 

Here P x (C 2 ), P 2 (C 2 ), P 3 (C 2 ) denote certain undetermined functions of U, V, in 
which these variables appear only in the form U 3 + V 2 + W 2 or G 2 . The first term of F 
is evidently of odd degree in U, V, W combined, and the second and third are of even 
degree; it is convenient to denote them by 0(ll, V, W) and E(U, V, W), when we 
wish to refer to the odd or even part of F separately.] 

It is easy to see that, in a uniform gas, f satisfies the necessary conditions (3), (4). 
In the non-uniform case these conditions require F to satisfy the equations 



(17) 



/ FdUdVdW = 0, 



(18) 



U/ F d\J d\l dW = V/ F dU dM dW = W/ F dU dW dYJ = 



Clearly the odd part (l), V, W) of F satisfies (17), and the even part E (U, V, W) 
satisfies (18), but not vice versa, so that these equations place certain restrictions on 
and E. 



§ 3. The Equation of Transfer of Molecular Properties. 

§ 3 (A) The rate of change of vQ, the aggregate value of Q (u, v, w) per unit volume, 
may be analysed into three parts, viz., that due to molecular encounters (which we 
denote by AQ), that due to the passage of molecules in or out of the volume-element 
considered, and that due to the action of the external forces. The equation 
expressing this analysis may readily be shown # to be 



(19) 



dt 



(,Q) = AQ- 2 



x,y,z 



— ( ~Ci) ~ x / dQ 

dx m" \du 



We may define AQ by the statement that ( AQ) dx dy dz dt is the change produced 
by molecular encounters during time dt in the sum 2Q taken over all the molecules 

in the volume-element dx dy dz : evidently 2Q = ^Q dx dy dz. 

If in (19) we make Q equal to unity, in which case AQ is clearly zero, the equation 
becomes 



dv 

dt 



dvUa . dvV 



+ 



, dvw 
+ ~3^ 



dx dy 

\dx dy dz 



dv , dv , dv 
dx oy - czi 



which is the equation of continuity. 

* Of. Jeans' c Dynamical Theory of Gases/ 
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Except under the differential sign we shall write u~\J>v = V,w = \N, since we are 
supposing that u = v = w = at the point (x, y, z). The last equation consequently 
reduces to 

(20). i| = _(^ + |io + a^ 

v ct \ax oy dz 

In taking mean values of functions of U, V, W, as in (5), we shall neglect the part 
F in the velocity-distribution function f, in cases where the mean value is to be 
differentiated or multiplied by a small factor, since the resulting error is only of the 
second order. # Thus, in such cases, we shall write 



(2.1) U 2 C 2 * = £C 2( ' +1) , U 4 ^^ - }G 2{s+1 \ VW^- 1 * - T ] 5G 2(S+1) , 



(22) G 2s = 1 . 3 . 5 . . . {2s + 1) (2hm)- 9 

while, if either p, q or r is odd, 



(23) U^V g W r = 0. 

Since the equation of transfer involves derivatives of the first order only, it is 
sufficient, whenever the mean value of a function of u, v, w is to be differentiated, to 
expand it by Taylor's theorem in terms of u , v 09 tv , so far as the first degree only ; 
if, then, the coefficient of u , v 09 w Q is of type (23), the corresponding term may be 
omitted altogether. 

Case I. Q = u (u 2 + v 2 + w 2 ) 8 . 

§3 (B) When Q = u(u 2 + v 2 + w 2 )*> according to the principles just laid down we 
have 

9 / -TTx 3 



(^Q) = ^- .(U 2 C 2s + 2^ UC 2s + 25U 2 KU + t; V + ^oW)C 2(s - 1) +...) 



dx ' dx 



i 3 / n^T)\ _ 1.3.5.,.(2g + 3) 3 / 1\ 
3 S^° ' 3 fa V \2hm) 

1.3.5 ...( 2* + 8) / J_V j _A_^ +v A ( _1_\ +sv A f_L 



2hm) {2hm dx dx \2hml dx \2hm/j ? 



#-(^Q)='o, f-M})==o, 

oy x dz 

* Except in gases of very low density. 
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neglecting in each case quantities of the second order. Since Q is of the first order, 
to the same degree of accuracy as in the above equations, we have 

Again, writing u=z\J,v = V,w = \N after differentiation, we have 

'^ = & + 2s[)*(?<'- i > = i(2« + 3)C 5 = ~- ' 5t 3 ^ 2S + ^ (^~) S ' 



Here we have omitted UV0 2( *" 1) and UWC 2(s-1) , since when multiplied by X, which 
is of the first order, the result is of the second order, and hence negligible. 
Similarly 

The equation of transfer consequently takes the form 

(24) AUC 2? = 1 >3. 5...(2s + 3) / 1 W 1 dv 8 / 1 \ v x d I 1 \ 

3 \2hmJ [2hm dx dx \2hmJ m dx \2hmJ 



When 5=0, this becomes 



1 dv t 9 / 1 \ * 



2 hm dx dx \2hmJ m 

Now mAU is the rate of change of momentum per unit volume due to the 
molecular encounters, and, since action and re-action are equal and opposite, this 
change is zero. Hence we have, remembering that (2h)" 1 = ET = p/v, 



(25) pX = 



2h dx ' v dx \2hJ dx \2h) dx ' 



which is one of the equations of pressure of the gas. 

On substituting the value of X given by (25) into (24), the latter may be written 

W 1.3.5...(2, + 3),; AUb ~Tto' 

where we have used the equation (cf. 8) 

9h d ( 1 \ _ 1 8T 
Zfim dx \2hml " T dx " 

There are two equations similar to (26) giving AVC 2s and AWC 2 * in terms of 3T/3y 
and dTfdz. 

VOL. CCXVI. — A. 2 E 
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Case II. Q = u 2 (u 2 + v 2 + w?) s . 
§ 3 (C) Making approximations and reductions as in Case I., we have 



3 

a7 



s n 



(,q) = ± a,c^») 



1 . 3.5...(2s + 3) 3 / 1 

= z — ' — v — 7- 



s+l 



ot \2hin, 



1 .3. 5...(2s + 3) / 1 \ ,+1 \ldv,, ,,\13T 



3 



[v ot 



T dt 



3 



3 



( vU Q) = ^„{U 3 C 2s + 3M U 2 C 2 ' ? + 2sU 3 ('« U + t' V + w W)C 2(s - 1) +...} 

0$/ O.X 



^ ^ (C 2(s + l) + l5C 2(s + l) ) = - * - * - — (^ g+ 5 ) ,. / 1 V +1 ^r 



„ _^- 



\2hmJ dx * 



h {vvQ) 



dy 

3 



i/ {U 2 VC 2s + t\,U 2 C 2s + 2sU 2 V KU +v V + w u \N) C 3( - 1) + ...} 



/ 



1 



fJ „n /ipb+i) i 2 „n2( s +D\ _ i . 3 . 5 ... (2s + 5) / 1 y +1 3y 



.d2 



1 , o . o . . . ( Zs -f- / / I 

" y 



15 



2Am, 



9w 

82; 



3Q N 



„^1 = 2UC 2s + 2sU 3 C 2(s - 1) = 0, 

ou I 



dv 



0, 



3Q 

dwj 



0. 



The equation of transfer may therefore be written 



(27) 



AU 2 C 



2rt2s 



1.3. 5. ..( 23 + 3) / 1 



\s+l 



15 



2hm, 



5 \v 3t + ^ +l ^T3* 

+ ( 2s + 5)(3^ + ^ + 3w 



3# dy 



dz 



When s = 0, this becomes 



AU 2 - 



j/ 



2 km 



13^.1 dT . n dn 



v dt T dt 



dx 



dv {) ?a.v 
-+- -+- 



% 



^ 



If to this be added two similar equations giving AV 2 and AW 2 , on the left-hand 
side we have A (U 2 + V 2 + W 2 ), which is the rate of change of molecular energy due to 
encounters ; by the principle of conservation of energy this is zero, so that 



~ , 1 dv 1 9T 

^'vdt T dt 






+ 5 (it 

\OX 







dy 



02 



0, 



or, by (20), 



(28) 



i^ar 
t a* 



2. 
3 



f 'dll dv 
dx dy ' dz J 



1 01/ 



CW \ _ 2 
-r — I — - 3 



v 



dt 



/9 06 
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On integration this gives T/o"" 8/s = constant, or, since p = Ri/T, it is equivalent to 
pp~* ls = constant ; this is the law of adiabatic expansion for a gas which possesses only 
translational energy. 

Eliminating - -- and = -—- from (27) by means of (20) and (28), we have 



(29) AU 2 C 2 * = l-3-5-(2* + 3) /l j' +1 f_«( 2 5 + 5) & + ^ + 



8u' {) 



15 \2hmJ 

\ ox ay oz I 



dx dy dz J 



— 1-3.5... (2s + 5) 9 / 1 V +1 (o 3^0 dv () 9W 
45 \2hmJ \ dx dy dz , 



or 



/o \ (2/?m)* 41 45 A . .apar = o d^o 3t\> thgp 

W/ 1.3. 5. ..(2^+5) 2v ox dy dz ' 

By transformation of axes, or otherwise, we m&y deduce the equation 

(31) - (f^: 1 -r f A2VWC- = 8 (^ + ^') . 

1 .3 . 5 ... (2,9 + 5) 2v \dz dy J 

§ 4. The Effect of Molecular Encounters. 

§ 4 (A) In this paper our primary concern is with simple gases containing molecules 
of one kind only ; the difficulties are much enhanced when molecules of two kinds are 
present, especially as regards the equations of transfer, and the final determination 
of the coefficients of F when AQ has been calculated. These matters will be dealt 
with in a future paper, on diffusion and the general theory of composite gases. In 
the calculation of AQ, however, there is scarcely any advantage in making the 
restriction to one kind of molecule only, and it is convenient to carry out the 
calculation for a composite gas in order that the results may be quoted without 
repetition in the later, more general, investigation. 

The notation of § 2 may be adapted to the case of a composite gas without further 
change than the addition of suffixes 1, 2 to denote to which group of molecules 
a symbol such as v, m, U, V, W, /, F refers. The mean velocity of the. two groups 
will be supposed the same, so that (u 0> v 09 iv ) is the same for both, either 
separately or together ; similarly, their temperature or mean energy, and their 
relative densities (vjvg) are supposed constant. All the remarks made concerning 
f and F hold both for f x and F 1? and f 2 and F 2 , these being functions respectively 
of (U l3 V 1? W^ and (U 2 , V 2 , W 3 ) ; they may each now be expected to involve v 1 :v ii 
and m l : m 2 in addition to the quantities mentioned in § 2. A further important 
consideration which did not arise there is that f x and f 2t or F x and F 2 are similar, 

Jj R <_/ 
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in the sense that either may be obtained from the other by interchanging the 
suffixes 1, 2. 

Notation for an Encounter. 



s- 



§ 4 (B) Before proceeding to the actual consideration of the dynamics of an encounter 
between two molecules m u m 2 , it is convenient to explain the notation to be used. 
The symbols m 09 /x 1? //<>, /x 12 , and /x 21 are defined as follows :— 

(32) m = w^ 1 + •r/^ 2 , 



(33) mi = Wi/m , fM 2 = m 2 /< 



mjm 



0? 



(34) /x ]2 = mi/ra 2 == ^//* 2 , M21 = mg/m-x = ^/mi, 
so that 

(35) M1 + M2 ™ 1? Mi 2 + M2 2 — 1 — 2//J/X2, 

Velocities will be represented either by their #, '?/, 2; components or in vector 
notation. The components of the actual velocities of the molecules will be written 
(U, V, W), while those of other velocities, such as the velocity of the mass-centre G, 
or the relative velocity, will be written (X, Y, Z). The amplitude of a velocity will 
be denoted by C, and the vector itself by the same symbol in small type with a bar 
beneath, viz., c. 

The velocities of the molecules m l9 m 2 , and of G . will be distinguished by the 
respective suffixes 1, 2, 0, while the suffix R, similarly, will indicate reference to the 
molecular velocities relative to G or to one another. As regards time, square brackets 
enclosing a symbol, such as [X ], [ft], will indicate reference to some particular 
(arbitrary) instant during the encounter ; a symbol without brackets but with an 
accent (') will refer to any instant after the encounter,* while when there is neither 
bracket nor accent it will refer to any instant before the encounter. 

Analysis of the Motions in an Encounter. 

§ 4 (C) In -the above notation the initial and final molecular velocities are respec- 
tively ft, c 2 and c\, c' 2 , or (U l5 V, W^, (U 2 , V 2 , W 8 ) and (U ; 1; V' 1? W\) 9 (U' 3 , V' 2 , W' 3 ) ; also 

(37) C 2 - U 2 + V 2 +W 2 , 

* That is, any instant after the molecules have separated beyond the distance (which in actual gases 
is, at most, very small) at which their inter-action is appreciable ; the words " before the encounter " are 
to be interpreted in a similar sense. In this sense the velocities of the molecules before and after the 
encounter are definite and constant. 
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where 0, U, V, W all have the same suffix 1 or 2, with or without an accent ('). 
Similarly, the mass-centre G has the velocity c or (X , Y , Z ), and since (by the 
principle of conservation of momentum) this remains invariable throughout the 
encounter,* we have 



(38) 



Qo = Mi£i + fi£'2 = fi\Q f \ + ^£ f 2 = Mi [&] + M 2 [&_ 



or 



Since, by (38), 



c 



c 



I/O — 1/ — L^OJ 



C r 



(39) m x {[cj-fco]} = -m 2 {[c 2 ]-[c ]} = m 0y a ]y u 2 {[&" 



L&J/ 



ra (miM 2 ) 1/2 [ftt" 



where [c R ] is defined by the equation 



(40) 



[C K ] = Gi 1 /* a )M[&]-[&]}> 



we see that the momentum of the molecules, relative to G, is equal in magnitude but 
opposite in direction in the two cases, its value being ±m (^t 1 ^ 2 ) 1/2 [ci]. The relative 
velocity of the two molecules is, by (40), equal to (m^)"^ [c k ] ; this varies throughout 
the encounter, owing: to the inter-action of the molecules ; its initial and final values 
are given by 



(41) 



Qr 



\/*l/*2/ \£l"""£2/j 



Sk 



U1M2) 1/3 (e'i — ^'2)3 



which are special cases of (40). 

Equations (38), (41), and the reciprocal equations 



(42) 



l k 



Q\ "~ &) + M21 £r 



V. 



^2 — £o~"~Mi2' ^lb 



(43) 



Q 1 — ±0 + M2.1 t' 11 L' 2 — i'O Ml2 ^ ]t? 



indicate that 0, c 2 or c' 1? c' 2 are equivalent to c , c K or c , e' E , as specifications of the 
initial or final velocities of the molecules. Hence the problem of determining the 
final velocities of two molecules after an encounter, in terms of the initial velocities 
and whatever further independent variables are necessary to define the encounter, is 
equivalent to the determination of c' R in terms of c R and the variables of the 
encounter. Thus, in consequence of the invariability of c , the velocity of the mass- 
centre, we need only consider the motion relative to G, i.e., the motion referred to 
uniformly moving axes with G as origin. 

* We here suppose that the effect of the external forces during the brief interval of encounter is 
negligible ; this is legitimate if the gas is " nearly perfect " (cf. § 2). 
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The Motion Relative to the Mass-centre. 

§ 4 (D) Relative to these axes the molecules are initially moving along parallel lines 
with equal and opposite momenta ±m ( y u 1 ju 2 ) 1/a c K , by (39). The plane containing these 
two lines is clearly the plane in which the inter-action and motion of the molecules 
will take place during the encounter. It is parallel to c E , but its orientation e about 
this direction is independent of c , c Ji9 i.e., it is one of the additional variables needed 
to specify the encounter, and, similarly, so also is the perpendicular distance jo between 
the initial lines of relative motion. It is convenient to measure e from the plane 
containing c and c u . 

In the plane of relative motion so defined, the molecules describe orbits which are 
similar to one another (the origin G being the centre of similitude), and symmetrical 
about the line of apses (i.e., points of closest approach). Each orbit has two 
asymptotes, one being the initial, the other the final line of motion ; the distance 
between the pair of final asymptotes is clearly equal to that, p, between the initial 
asymptotes. The angle Xi2 between the two asymptotes of either orbit measures the 
deflection of the relative motion due to the encounter ; for molecules of given types 
it is a function of p and C K # only, the nature of the function depending on the law of 
inter-action between a molecule m 1 and a molecule m 2 . We shall find it convenient, 
for the sake of generality as well as of brevity, to retain X12 as an unspecified function 
of p and C K in our equations ; the special properties of the molecules under 
consideration are, throughout our work, involved only through the dependence of 
X12 on p and C R . 

It is easy to see that the magnitude of the relative velocity (mi^)"^^ * s unaltered 
by the encounter, i.e., 

(44) C a = C' B : 

for by the equation of energy we have 

(45) i(mA s +™A 2 ) - iKcv+^c/) = W<V+(V) - W(V+cy) 

by (42) and (43). 

The Velocities in Spherical Polar Co-ordinates. 

§ 4 (E) The above analysis of a molecular encounter may be made clearer by the 
following figure, in which x, y, z, c , c K , c' K are the points in which a unit sphere 

* That is, on p and on the amplitude C R of the vector c R . 
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centred at O is cut by radii parallel respectively to the co-ordinate axes and to c , c H , 
and c' R . Then we have 



(45a) 



€ Gnt/i?C 



U u )l° R? 



Xl2 



GrOc R* 




It is convenient also to use certain spherical polar co-ordinates, as follows, taking 
Ox, Oxy as initial line and plane for c R and c' R , and Oc R , Oc n x 9 or Oc' R , Oc' R aj for c . 
Thus we write 

(46) R = CflOa, fl' R = c' R Oaj, = c Oc R , 0' o = tf Oc' R , X = c Ox, 



(47) 



0R — C R X i/9 R ~ C ll X U^ 00 — C C R^> — C C K^« 



Evidently we have 



(48) 



cos 0' o = cos O cos xis+ sin O sin x i2 cos e, 



(49) 
(50) 



cos 6' R = cos 6 n cos X12+ sin 6 n sin xi 2 cos e+ O , 



cos X = cos O cos R + sin O sin R cos <£ 



0> 



= cos 0' o cos 8- n + sin 0' o sin 7 R cos <j>' . 



Expressions for the Velocities After an Encounter. 

§ 4 (F) We have thus indicated how the final molecular velocities c\, c' 2 are to be 
determined (cf 43) in terms of the initial velocities c u c 2 or c 0i c R together withjp and € 
(these being the eight independent variables of an encounter). This has been done 
by showing how c' R depends upon c R , p and e ; it has in fact been shown that the 
spherical polar co-ordinates of c 7 R , referred to c R and the plane e , c R as initial line and 
plane, are C R , X12 ( a function of p and C R ) and e. Hence we may at once write down 
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the expressions for the initial and final velocities and velocity components in terms 
of c , c H , p and e, as follows :~ - 

(5i) mA 2 = e 12 (cos O ), m 2 c 2 2 = e 21 (cos e ), MiC? = e 12 (cos e' ), ^c? = e 21 (cos e' ), 



(52) < 



ya/^Ui = juli^Cq COS X + H^Cr COS k , M2 i/: 'U 2 = M 2 1/2 C COS X — /ULi^G j> cos K , 

Mi 1/a U'i = mi 1/s C cos \ + jut.2 l3 C R cos 0' K , m 2 1/2 UV= M2 l/i C cos X— ^i^Cp, cos 0' K , 



where we have adopted the convenient notation^ defined by 

fe r ,(cOS 0) = /z 1 C 2 + /x 9 C K 2 +2 ( l a 1 ^ 9 y k G {} C R COS 0, 

(53) 4 

[G 21 (cos 0) = yu 2 C 2 + yuiC K 2 — 2 (m.iM2) 1/2 C C k cos 0. 

Equations (51) to (53) are merely particular cases of (42), (43), expressed in terms 
of amplitudes (51, 53) and of x-components (52). The latter might also have been 
written in terms of the components of C and G Ry as, for example, 

(54) \J\ = X W /j X' R = X W 9 {X B cos xia + (Y tt 2 +Z R 2 ) l/3 sin Xw cos (e + ^ )}, 

by (49), writing (X , Y , Z ), (X R , Y K , Z R ), (X' R , Y' H , Z n '), for the components of c , c K r/ K . 
Equations similar to (53), (54) may easily be written down also for the y and z 
components of the velocities. 

The Dependence of U'i, V' 3 , WVon xn> 

§ 4 (G) From (51) and (54) it is clear that any function Q 1 (U'i, V'i, W'i) of U'i, V\, W ; 
is a function of U x , V 2 , W x , U 2 , V 2 , W 2 , jp and <?, or of U 1? V 1? W 1? U 2 , V 2 , W 2 , X u an cl e, 
since p is involved only through X12 (though X ia * s not entirely independent of the 
preceding six variables, since it depends upon G R ). If Qi (ll'i, V'i, W'i) be regarded as 
a function of xis, when xu is- made equal to zero it reduces to Q x (Ui, Vi, Wi) simply : 
this may be seen either from (5l)~(54) or, still more readily, from the figure on p. 293, 
since when xu = 0, c\ t becomes identical with c m and hence by (42), (43), so also does 
e\ with c a . 

Transformation of Co-ordinates. 
§ 4 (H) In § 5 we require the Jacobian of transformation 

j a^v^w 7 !, u^v^w',) 

3(Ui, Vi,W 1? U 2 ,V 2 , W 2 ) 

between the initial and final velocity components, p and e being constant. Since the 
motion during an encounter is reversible, the relation between the two sets of velocity 

* In § 7, for the sake of brevity, we shall write 0i 2 , 2 i, 6' 12 , 0' 2 i respectively for 12 (cos ), 2i (cos # )> 
012 (cos 0' o ), and 2 i (cos #' ). 
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components is reciprocal, so that J = ± 1. It may readily be seen that the positive sign 
is the correct one, by considering a particular case of variation, say dl)i = d\J 2 = d(J, 
dVi = dV 2 = c?V, dSNi = c£W 2 = c?W. This is equivalent to the addition of a small 
velocity (dU, d\l, cA/V) to the whole system; obviously this will reappear in the final 
velocities, so that also d\j\ = dl)' 2 — dl), dV\ = c£v 7 2 = d!V, dW'i = c?W 7 2 = d\N. In 
this case, and therefore always, we have J = 1, so that 

(55) dll'i cN\ d\N\ d\i f 2 dV f 2 dW 2 = c?U x dVi d\N 1 dU 2 dV 2 d\N 2 . 

With a little more trouble this might also be proved analytically from the equations 
of this section. 

From the component equations corresponding to (42), (43), i.e., from 



(56) < 



Ui - X + M21 1/2 X E , V x = Y W'Y* Wi - Z + /x 2 ; /2 Z 

U 2 = X — ycc 12 2 X R , V 2 = Y o~~"Mi2 ! 'e> W 2 = Z — M12 2 Z 



2 
11 5 



(57) C 2 = X 2 + Y 2 + Z 2 , C B a = X s 3 +Y B 2 + 

it is easy to prove that 

(58) aOji^w. u 2! v,, w.) = _ M _> K 

v (Xo? Y , Z 0? X r , Yr, Z r ) 
Hence, by further transformation to polar co-ordinates, we have 

(59) c?Ui d\/ 1 d\Nx d\J 2 d\l 2 d\N 2 = - (miM 2 )" 3/2 ^X c?Y c£Z c£X r c?Y e dZ K 

= — (mim 2 )~ 3/2 C 2 C r 2 c?C c£C r c£ cos # c? cos 6 n dcp d(p n . 

Since dUid\/ 1 d\Nid(J 2 dy 2 d\N 2 is essentially positive, the negative sign on the right 
of (59) must be made positive, if the limits of C, cos 0, and <j> in each case are taken 
as to + o°, — 1 to +1, and to 2ir respectively; it may readily be seen that the 
negative sign corresponds to reversed limits of integration of one of the variables cos 9. 

§5. The General Expression for AQ x . 
Definition of AnQi and A 12 Q lr 

§ 5 (A) The rate of change of ^Qi due to molecular encounters, i.e., AQi, may be 
divided into the two parts AnQi, A 12 Q X due respectively to the encounters of the 
molecules m 1 among themselves, and those with molecules m 2 . Thus 

(60) AQ 1 = AnQi + AjaQ^ 

We shall chiefly consider A 12 Q 1? whence AnQx may be obtained by changing the suffix 
2 into 1 throughout. 

VOL. OOXVI. — A, 2 S 
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The Expression for A 12 Qi. 

§ 5 (B) The number of molecules m ± having velocity components lying between the 
limits (U l5 Vi, Wi) and (Ui + dUi, V 1 + d\f u Wi + dWj) is, by our definition of f x (U 1? Vi, Wj), 
equal to 

^i/i(Ui,Vi, Wi)dUidVidWi 

per unit volume. The number of encounters in time dt of any one of these, with a 
molecule m 2 having velocity components lying between the limits (ll 2 , V 2 , W 2 ) and 
(U 2 + c£U 2 , V 2 +<iV 2 , W 2 + <iW 2 ), the variables _p, e of the encounter lying between p and 
p + dp, e and e + de, is equal to the number of such molecules m 2 contained within a 
small cylinder of length (miM2)~ 1/2 C r dt and of sectional area p dp de> i.e.> to 

V2 (miM2)" 1/2 / 2 (U 2 > V 2 , W 2 ) C r jd c£p de dU 2 dV 2 d\N 2 dt. 

Thus the total number of encounters of the above type, per unit volume per unit 
time, is 

(61) viv 2 (fx^y Vi (Ux, Vi, W^/a (U 2 , V 2 , W 2 ) C R ^ d_p ck cM dVx dWi c?U2 ^V 2 dW 2 . 
At each such encounter the change in the value of Q 1 (\J U Vi, Wi) is clearly 

(62) Q, (U',, V'x, W'J-Qx (U 1? V 1? WO, 

or Q 7 ! — Qi, as we shall write it for brevity. 

We shall include the effect of all possible encounters per unit volume per unit time 
if we integrate the product of (61) and (62) over all values of e (0 to 2tt), p (0 to go) 
and (Ui, Vi, W^, (U 2 , V 2 , W 2 ) (each from — oo to + oo). Such an integration will 
include encounters which are not binary, but our postulate that the gas is nearly 
perfect (§ 2) implies that our integral would be altered only inappreciably if the 
upper limit of integration for p were not infinity but equal to the very small distance 
at which two molecules cease to exercise any appreciable inter-action. Hence, 
throughout this paper, where no limits of integration are specified, it is to be under- 
stood that they have the above values. Thus we have 



r 



(63) A u Qi = v lV2 (mim 2 )-' /2 (Q'i-Qi)/i/ 8 C B i> dp de dUx dWi dWi d\J 2 dW 2 dW 2 . 



%i *> 



The termfif 2 in the integrand may be written 
(64) ^J' , ^y /, e-M- 1 o,. +wrftf , {i + Fl (U x , Vx, Wx) + F 2 (U 2 , V 2 , W 2 )} 

= W 2 ) 3/2 (^) 3 e- to » (c » 2+ V (l + F 1 + F 2 ), 

where, in the first line, we have neglected F^, which is a second-order quantity, 
while in the second line we have made use of (45). 
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The unit term in (l+Fi + F 2 ) may be omitted. 

§ 5 (C) It is easy to show that the part of (63) which arises from the unit term of 
(l+Fi + Fg) in (64) is zero. For it may be written in the form 



where 
</> = 



%s %s %s J %) 



<pviv 2f x lf i 2 (hm /ir) 3 C R _p dp de, 



{Qi (U'i, Vi, W'O-Qx (U 1; Vi, Wi)} e-'^^+V dUx dVj d\N x d(J 2 dV 2 c?W. 



Now by (45) and (55) the latter is equal to 



(65) 



v %> 



Qi (U'i, V'i, W'i) e - hlmifflt+m ^ A dU'i dN\ dW'i dU' 2 dv' 2 dW. 



V */ <J »/ «. 



Qi (Ui, Vi, Wi) e- fc <™i°i s +™*w cZUi dVi c?Wi c?U 2 c?V 2 dW 2 . 



But the latter two integrals are equal, since they are definite integrals differing only 
in the symbols used to denote the variables. Hence (65) is zero, and the unit term 
in (l-f Fi-j-F 2 ) may be omitted from A 12 Qi. 

The same result can be seen also in another way : the part of A 12 Qi under 
consideration is that obtained by putting F x = F 2 = mff 2l i.e., it is equal to the 
value of A 12 Qi in a uniform gas. In a uniform gas, however, as we may see from the 
general equation of transfer (19), A u Qi = A 12 Qi = 0, whence the result follows at 



once. 



If Q (U, V, W) is of odd degree, the even part of F(U, V, W) contributes 

nothing to A 12 Q, and vice versa. 

§ 5 (D) We may now, therefore, write A 12 Qi in the following form, transforming 
the variables (U 1? W u Wi), (U 2 , V 2 , W 2 ) to (X , Y , Z ), (xr, Y r? Z e ), by (56), (58). 

(66) A 12 Qi = ^ 2 /xiM2i — - ] 



\ 7T 



*/*}*. 



W2 ( 



MlMl) 



- l k 



(Q^-Qi)e-^ (Co2+c * 2) (F^ 

' hm Y 

(Q'l-QO^^^+^^Fx + Fa) C n p dp dedX dY dZ dX R dY n dZ R . 



We here suppose the functions Q (U, V, W) and F (u, V, W) expressed in terms of the 
new variables and (in the case of Q\) of e and X12 ( or p)> We are concerned both as 
regards Q and F only with terms which are integral in the variables U, V, W ; in 
reckoning their degree we shall make no distinction between Ui and U 2 , &c, or 

2 s 2 
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between X and X R , # &c. Then since the equations of transformation (56) are linear, 
any term U^V^W^U/V/W/ in Qi(Fi + F 3 ) transforms into the sum of a number of 
terms X a Y & Z %/Y R e Z/ such that 

l-Vp — a + d, m-\-q — h + e, n + r = c+jt. 

This is not true in the case of Q\ (F x + F 2 ), since (by 54) U'i, V ; i, W'i are not rational 
functions of the variables X, Y, Z, but it is true of (Fi + F 2 ) jQ'i de, since the integration 
with respect to e causes all the irrational terms in Q\ to disappear."!" This may be 
proved quite generally, but it will be sufficient here to indicate the proof for the case 
Qi = UiCi 2s , s being any positive integer. We may write 



Q\ = U'iC? 5 = (x + aC B .sin 6 R cos e + fo) (C 2 + 2aC C R sin 6 cos e) 5 , 
where 

a = fi 2 i k sin Xm X = X + [i 2 i%i cos X i 3 , 

C 2 = .C 2 + M2 A 2 + W' 8 (X X R + Y Y R + Z Z E ) cos Xl2 , 

so that X is of the first degree in X or X R , and C 2 is of even degree in the variables 
(X 0? X E ), (Y , Y r ), (Z 0) Z e ). The only terms in Q\ which do not vanish on integration 
with respect to e are of the form 

X { s G 2p (C 2 ) 8 - 2 ? (2aC C R sin 6 cos e) 2 ?} 
or 

(C C E 2 sin 6 sin 6 R cos 2 ^ +1 e cos H^ ) {2a 2 s G 2p+1 (C 2 )-*- 1 (2aC C R sin O ) 2 ^}. 
Now we have 

(C C U sin e o f = C 2 C K 2 (l-cos 8 O ) = {G 2 C U 2 -(X X R + Y Y K + Z Z B ) 2 }, 

which is an even function of X, Y, Z, and can be included under the symbol C 2 . 
Thus, on integration with respect to e, the above expressions become (apart from 
a factor not involving X, Y, Z explicitly) 

XC 2s , (C C R 2 sin 6 sin e R cos <p ) G 2{s - 1] 

and by (50) the latter may be written 

C C R 2 (cos A-cos 6 cos 6 R ) C 2( *- 1} - [XoO B 2 -X B (X X R + YoY R +Z Z E )] C 2(s - 1} . 

Both these expressions, and consequently JQ X (ll'i, V'i, W'i) c£e as a whole, are of the 
form XC 2s in the sense above defined. Similarly it may be shown that JuVCm 2 * de is 
even in all three variables (X 05 X B ), (Y , Y e ), (Z , Z r ). 

* So that, for instance, x 2 , x x R and x E 2 will all be regarded as even functions of x. 
f The explicit occurrence of x, y, z in f Q'i de is here referred to ; the latter may involve C R irrationally 
through xi2- 
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In the integrand of (66), the exponential term and C R (whether occurring as an 
explicit factor or implicitly in xu) &re even functions of X, Y, Z. Hence a term such 
as X a Y 6 Z c X/Y/Z/ in j(Q'i-Qi) (Fi + F 2 ) de will contribute nothing to A 12 Qi unless 
a, b, c, d, e,f and a fortiori a + d, b + e, c+f are all separately even. In view of 
what has been proved above, therefore, it appears that in Q (U, V, W) F(U, V, W) 
only the terms which are even in U, V, W separately contribute anything to A 12 Qi. 
Hence if Q is odd in U, only the part of F which is likewise odd in U need be 
considered, while if it is even in U, only the even part of F need be considered. 

Introduction of l x (xi 2 ). 

§ 5 (E) We now make the final transformation of A 12 Qi by adopting polar co-ordi- 
nates in place of (X , Y , Z ), (X R , Y E , Z R ), as follows :— 

(67) A 12 Q X = VlV2 {^ 2 y- (^f |||e"^(Co^ V ) {Ij ^^ (o)} q,^ dpdCo dCjb 
where 

(68) IxGte) = QiCU'lV'lW',) {F 1 (U 1 ,V 1 ,W 1 ) + F 8 (U 8 ,V a ,W 8 )}ckdcosMcosM0o<fyB,- 



J J J 



(69) 1^0) = Q, (Uj, Vi, WO {^(U^VlWO + F^Us,, y 2 ,VJ 2 )}dedGose dGosd li di> d<p Ji . 



J *> 



Evidently (cf § 4 (G)) the latter is obtained when xu is made zero in Ii (xi 2 ), since X12 
is not concerned in the integrations of (68), (69), being a function of p and C R only, 
while when X i2 = we have Qi (U'i, Vi, W'i) = Qi(Ui, Vi, Wi). Hence, in calculating 
A 12 Qi we shall concern ourselves only with Ii (x i2 ) until we come to the integration 
with respect top, C R , C . In so doing we shall, from the outset, omit from F (U, V, W) 
those parts which, in accordance with '§ 5 (D), contribute nothing to the final result. 

§6. The Form of the Function F(U, V; W). 

The two special forms of Q x which we consider are U^C^ 8 and (Jfi^ 8 ; the only- 
parts of F(U, V, W) which are relevant in these cases are respectively the part of 
Ej -f E 2 , which is even in V and W 2 , and Ox + 2 ; the notation here used is that of 
§ 2 (E), p. 283. From (26) and (30) we see that AUfij 28 involves the space derivatives 

of mean properties of the gas only in the form — — , while AU/Ci 28 similarly involves 

only 2 ^ - ^ - ^. We deduce from this that (U, V, W) must certainly include 

J dx dy dz x 

the term 
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and that E((J, V, W) must similarly include the term 

(71) (c 11 U 2 + c 22 V 2 + c 3 3W 2 + c 23 VW + C3 1 WV + c 12 UV)P(C 2 ), 



where 



r 



(72) 



< 



= 2 du dv {) dw ^ ^ == 3 fdv, + dwA 

dx dy dz 23 \ cz dy j ' 

3^0 dw t) du ___ Q fdw . 3^ ( 



22 "* 3?/ cfe dx ' 31 ' \ 3x 



z 



V 



j = o 9 ^° 8 ^° 9 ^° c =3 f ^ + ™° 

3z dx dy ' 12 \;3^/ 3oj/ * 



The factor of P(C 2 ) in (71) is equal to 3S / — C 2 S, by (12) and (15), and is therefore 
an invariant with respect to an orthogonal transformation of axes. 

Further since, by (26) and (30), no other derivatives of T and (u , v 0i w ) occur in 
AUxCi 28 or AU/Cj 2 ', we conclude that none such appear in F(U, V, W) — at any rate, 
to our degree of approximation; thus the other terms in (12)-(15), while they 
possess the invariant property, do not satisfy the other conditions which must be 
fulfilled by F (U, V, W). 

We therefore conclude that F(U, V, W) is composed only of (70) and (71) to our 
order of accuracy, and we shall suppose that the two functions P (C 2 ) are expansible 
as power series in C 2 . Throughout this paper we shall assume that all convergency 
conditions necessary for the validity of our analysis are satisfied ; the justification 
of this assumption would offer serious difficulty, and the investigation would lead 
us into regions of pure mathematics which are largely unexplored, and would be 
unsuitable in the present paper. In § 10 we shall see that numerical approximations 
for the most important molecular models confirm the assumption of convergence 
sufficiently for our purpose. 

It is convenient to write our expression for F (U, V, W) in the form 

(73) F(U V W) = -B ft -/u — +V— + W — ) 2 (2hm) r ^ ^ 2r 

\(6) * {u, v, vv; r> 0T ^u ^h-v ^ + vy ^j^ ^ 3 ^ 5 _ ( 2r + 3)r ^ 

-C 2^m(c u U 2 +c 22 V 2 + C3 3 W 2 + c 23 VW + c 31 WU+c 12 UV) 2 - K f L _ ± _^ y r G 2r . 

In the first line, when r = 0, the factor r in the denominator is to be omitted. 
The suffix 1 or 2 must be added to m, U 5 V, W, C, /3, y when we wish to distinguish 
between F^Ui, Vi, W,) and F 2 (U 2 , V a , W a ). 

Since, by (72), 

(74) c u + 033 + 038 = 0, 
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it is easy to see that (18) is satisfied by this form of F, while in order to satisfy (17) 
we must have 



00 



(75) /8_ 1 + 2 j 8V(r+l) = 0. 



The products B /3 r , B y r are quite definite, but B and O can evidently be assigned 
arbitrarily ; we shall decide that their values, though unspecified for the present, are 
alike for ¥ x and F 2 . 

The above expression for F (U, V, W) is equivalent to that obtained by Enskog (§ l), 
by an entirely different method, But the chief difficulty of our problem, and one 
hitherto unsolved, lies in the determination of the coefficients /3 and y ; this is effected 
in the present paper by means of AQ. 

§7. The Calculation of AQi. 

§ 7 (A) In calculating AQi we shall deal chiefly with A 13 Qi (cf. § 5 (A) and (67)). 
The particular forms of Q x which we shall consider are 

(76) Q, = (2Am 1 ) s+4 U 1 C 1 2s = 3&i (s> , 

(77) Qi = (2^m l ) s+1 U 1 2 C 1 2s = d w . 

In accordance with §5 (D), the only part of F(U, V, W) which is relevant to 

ABi (s) is 

/itq\ t> 1 8T . .^ (2Am) r . r. ™ r 

(78) _J *°T3b o 1.3.5. ..(2r + 3)r*- lU ' 

while that which alone concerns Ad (s) is 

(79) -2hmC (c u U 2 + C22 V 2 +C33W 3 ) I Y^f^J^E) ^ 

As to the latter, since the remainder of the integrand of ACi (s) is symmetrical with 
respect to V and W, the parts of this integral arising from V 2 and W 2 in (79) are 
equal, so that c 22 V 2 + c 33 W 2 can be replaced by i(c 22 + c 33 ) (V 2 + W 2 ) = — Jc n (0 2 — U 2 ), 
by (74). Hence for our purpose (79) is equivalent to 

(so) -*( 2Aw )°» c ^ 3u2 - C2 )|r7^r|^^^ c " 

We shall denote by &i 2 (nsi) the part of A 12 Bi (s) which arises from the term 
— (2^m 1 ) r+i UiCi 2r in Fi(Ux, Vi, Wi), and by 6 12 (?Vi) the part arising from the corre- 
sponding term of F 2 (U 2 , V 2> W 2 ), in each case the numerical and other factors in F 
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being excluded.* The corresponding portions of Ii (X12) will be denoted by I (n^i, X12) 
and I (r^Sx, X12) respectively. Then 



(81) I(ria 1} xia) 



*> %j 



(2A / u 1 m ) r+s+1 U / iC? s UiCi 2 '' de d cos 9 d cos B d«£ dfa 



(82) I(r s s 1)Xl3 ) 



(2/im ) r+i+ V + W + * U / iC 7 1 2s U 2 C 2 2r de d cos O d cos B d«£ d«£ 



Bi 



J J J V 



(83) 6 ia (r lSl ) = - ^ ( MlMa )-v.fey 



7T / 



** J 



e- tao(w+V) {I (rx«x, Xl2 )-I (n«x, 0)} C 8 C B 3 JP # ^C dC B , 



(84) 6 u (r A ) = -^.(w)- ,fc ^Y 



\ 7T / 



J J 



e -WW + «y> {I (rA> ^.j (rA> 0)} Co3(Vp # rfCo ^ 



The similar quantities relating to A^®^ will be denoted by c 12 (r 1 Si)> c 12 (r^) and 
J ( r i s i) Xm)> J (?Vi, X12) respectively, so that 

(85) J^, X12) 



«/*'•. 



2 



»/ %/ 



(2V 1 r« ) r+ ' +a U?C? S (3U/-C! 2 ) C * de d cos O d cos B d^» cty K) 



(86) J(r 2 Si, X12) 



*/ «v */ 



1 (2^)'+'+ V + V/ +, U7C7' (3U/-C/) C 2 2r de d cos O d cos B d^> cty H , 



(87) Cj^r^) = -^ 



_i /2 /Am \ 3 



"l"2 VM1M2) V ' 2 ( J 



J J J 



e-*- w+0 -*> {J (r A> xt 2 ) -J (r A , 0)} C 2 C B 3 p dp dO dC B , 



(88) c 12 (r 2S] ) = _^ 3 (^ 3 )-'/,feY 



\ 7T / 



e - ftw .«v + o^ { j (rA) X]2 ) _ j {rsSu 0)} 0^3 ^ ^ dG ^ 



* We have here included a factor {2hm) 1 i% which does not occur in F ; this will be allowed for 
subsequently. 
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It 

The Integration ivith respect to 9 lb H , <p . 

§7 (B) In I(rs 9 X12) and J (rs, X12) it is clear from (48) to (53) that <f> n does not 
appear at all in the integrand, while 9 n and c/> occur only in the products U'iU and 
U'i 2 (3U 2 — C 2 ) respectively. We have 

U'iUi = G 2 cos 2 X + /x 2 ](V cos 9 R cos 6\ i + At 31 1/s C C K cos X (cos 9 n + cos 0' K ), 

\J\U2 — C 2 cos 2 X— C K 2 cos K cos O'r + CqC^ cos X (/x 2] 1/2 cos 0^ — Mi^ cos E ), 

and, remembering the values of cos X and cos 0' R , we have 



U'jUi d cos K d<f> dfa = §7r 2 {C 2 + /* 3 iC B 2 cos xi2~M2i 1/2 C C K (cos O + cos & )}, 
\J\U2dc0B Qjidfodfa = -|tt 3 {C 2 — C e 2 cos xi2 + C C R (/x 21 1/2 cos 0' o — /x 12 3/a cos O )}. 



»/ j 



j j 



In the notation of (53) the latter two equations may be conveniently re- written as 
follows : — 



/* s% 



Mi 



U'iUj cl cos 9 u d^dcj> ll = |tt 2 {0 ]2 (cos o ) +0 ]2 (cos 0' o ) -2^0/(1- cos X12)} 



(/Xj/X 2 j 



%/ & 



U'lUa d cos R d0 o c^ = |tt 2 W /2 02i (cos o ) +/x 31 1/s ia (cos 0' o ) 

+ 2 (/xiM2) 1/2 C K 2 (1 - cos x ) - (miM2)~ 1/2 C k 2 } 



Substituting in (81) and (82), we thus have 
(89) I(r A , X12) = $** (2hm Q ) r + t+l ff {e 12 +9^-2 M2 Cy (l~ cos Xl2 )} O'^O^dedcoBd, 



(90) I (r 2 .s 1? X12) - l^r 2 (2/m* )' + * +1 j J U\ + /iX + 2 (aW'<V (l - cos X]a ) 

— (miM2)~ 1/2 C k 2 } O / 12 *0 12 r de.d cos O . 
§ 7 (C) In the case of J (rs, X12X we have 

(91) ^UVUx 2 = (Mi^ cosXW<3*cose^ 

in which (c/ the figure on p. 293) X = c Ox, R = c K Oas, 0' B = c' R Cte. In the integration 

over the sphere, with respect to H and <£ , since O , e, X12 are constant the triangle 

CoCrc'h preserves its form, so that we may, if we please, regard x as the variable 

point and c^c'^ as fixed. Now it may readily be proved, by the method of 

"poles" in the theory of harmonic functions, that if A, B, C are three fixed points 

on a unit sphere, and P a variable point, then the integral over the spherical 

surface of 

cos 2 PA cos PB cos PC 
is 

t Vtt (2 cos AB cos AC+ cos BC). 

VOL. COX VI. A. 2 T 
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Applying this result to (91), identifying A, B, C with one or more of the points 
c o> c k> ^'li? an d P with x, we may with but little difficulty prove that 

*** \\\ U?U * d C ° S dn d(t>0 d ** = T%w2 { ° r * + 4Gl29/l2 + 9 '» ~~ 4m2 ° e2 ( ° 12 + & ^ ( l ~~ C ° S Xls) 

+ WC B 4 (l-cosxi 2 ) 2 }. 



Similarly we may show that 



Mi 



so that 



{%/%/* 



j j 



U? d cos R d$ d<p K = §7r 2 e\ 2 , 



Ml 



2 



u? (3u, a -cv) ^ cos B cfy ety E = 4-Tr 2 {e 12 2 +-fe 12 e' 12 +ey 

- 4/* a C M 8 (0 18 + G'u) ( 1 - cos xw) + WC K * ( 1 - COS xi 3 ) s } • 



Hence we have 
(92) J (r lSl , Xl2 ) = i^ (2hm Y^ 



{e 12 2 +fe 12 e' 13 +e' 1 /-4 M2 c E 2 (e 12 +e / 12 )(i- cos X]2 ) 

+ W(V (1 - cos X i 2 ) 2 } 0' 12 s 9i/ de tZ COS O , 
and it may be proved in a similar manner that 



./ 2 

i 2 



(93) j ( % Xl2 ) = iv (2/m )" +s+2 J j [ Ml2 e 21 2 +§e 21 e' 13 + / * itt e' 1 . 

+ 2C B a (//. 12 1/2 e 21 W*9'i 8 ) {2 (m,m 2 ) ,/2 (1 - cos Xl2 ) - (m^ 2 )- 1/s } 

+ cy {2 M 1 '* (1 - cos Xl2 ) - (^a)- 1 ' 2 } 2 ] e' ]a »e 21 ' & d cos o . 

The Expansion of(p 2 + <r 2 —2po- cos 0)" in a Legendre's Series. 

§ 7 (D) In order to effect the integration of I and J with respect to e and O we must 
have recourse to the expansion of 



(94) 



P„(/3, <r, COS 9) = (p 2 + cr 2 — 2p<r COS 0)" 



in a series of Legendre's functions. In a recent paper* I have shown that 



(95) 



n 



P n {p, <r, COS 0) = 2 (- 1)* (2£ + 1) B A /£ (p 2 , cr 2 ) P A .(coS 0), 



where P A . (cos 0) is the ordinary Legendre's function of cos 0, of type k, andf 



(96) 



"A* (p\ a 2 ) 



- P 



,k 11 



2 



Jh \^ + 2/t-k 2{r-t) 2i 

P (J 



k n 



(T 



2 



n* 



p! t=k (t + 3s) t (t — k) 



\ n + 2/t-k _2(r-t) 2t 

p • 



(T 



t-k 



* Chapman, ' Quarterly Journal of Mathematics,' p. 16, 1916. The expansion is there not limited to 
integral values of n, though these are alone considered in the present paper, 
t The constant k is necessarily a positive integer ; if k > n, n A h = 0. 
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In the last equation the symbol^, where q is integral, is defined thus : — ■ 

(97) p q = p(p-l) (p-2) ... (p-q + l). 
From (53) it is clear that 

n 

(98) {6 12 (cos e)} n = P B ( Ml C 2 , mAo - cos 6>) = 2 (2*+ 1) 'A 7 < 12 P* (cos 0), 

(99) {e a (coB d)} n = P„ ( M2 C 2 , miC k 2 , cos 0) = 2 (~l)*(2£ + l)'A* 21 P*(cos 0), 

ft=0 

where we have written, for brevity, 

(100) »A* 12 = "A" ( Ml C 2 , M2 C K 2 ) »A* 21 = -A* (mA 2 , M1 C K 2 ). * 

In our expressions for AQi, 9 takes the values 6 and 0' o , and the variable e 
is involved only through the latter angle, which occurs in Q\ 2 or 9 12 (cos 0' o ). 
In the expansion of the latter (cf. 98) in terms of P A (cos0'o), or, by (48), of 
T? k (cos6 cos xi2 + sin o sinxi2cos e), we shall make use of the following well-known 
formula in the theory of spherical harmonic functions : — ■ 

(101) P /c (cos 0' o ) = P A (cos 6 cos xi2 + sin 6 sin X12 cos (?) 

- p ; ;(cos e ) p a . (cos X ] 2 )+2 2 W— yf; p* 1 (cos e) p* z (cos X]2 ) cos u. 

TAe Integration with respect to e and . 

§ 7 (E) Since the integral of cos le with respect to e, between and 2tt, is zero unless 
Z=0, from (98) and (101) we deduce the result, 



C2tt n 



(102) \Q' l2 n de = 2wk (2k + l) w A\ 3 P fc (cos O ) P* (cos X12). 



Now from (89), (90), and (92), (93), it is evident that as far as concerns integration 
with respect to e and 9 we have to consider a number of terms such as 



(103) 



G m G\ 9 n de d cos d, 



12 w/t U/ V,WO l/Q j 



where 6 m may have the suffix 12 or 21, while. 0' n always has the suffix 12. Now Q m 
does not involve e, so that (102) suffices for the integration with respect to e, and 
leaves us with 



•1 [ m 

1 2 

-1 l/t=0 



(104) 2tt 2(±l) ;s (2£ + l) m A A 'P,(cos0 o ) 



1 2 (2*+ 1) "A*, S P* (cos 6 ) P* (cos X12) f d cos O 

2 T 2 
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in place of (103) ; in the first bracket the ambiguous sign is to be + in the case of 
12 w , and — in that of 21 w \ 

By the theory of Legendre's functions we have 



i 



(2/c+l) P A (cos d ) P z (cos O ) d cos 6 = if h ^ I. 



(2i + l) {P /c (cos0 o )} 2 dcos0 o = 2. 
J_i 



Consequently 
(105) 



m, n 



e^G'w* de d cos 6 = 4tt 2 (2k + 1) m A* 12 "A* 12 P 4 (cos X)2 ), 



i=0 



(106) 



?», n 



e 21 m e\/ de d cos 6, = 4t2 (-1)*(2A + 1) mi A* 21 b A* 12 P a .(cos X i 2 ), 



fc=0 



where the upper limit of h is the lesser of the two integers m, n. 
Applying these results to (89), (90), (92), (93), we have, therefore, 



(107) I(r A , X i a ) 



1 6 nr-3 

- ; v-7T 



!tfft 



>o) r+ * +1 



r+1, s + l 

V J" ^+ 1 A k s k k l r A /c «+l A /c 
Z/ t ^ 12 -^ ]2~r ii- 12 -&- P 
k = 



12 



2 M2 CV A 12 S A* 13 (1 -cos Xia )} P, (cos Xl3 ) ; 



(108) I (r 2Sl , Xl2 ) = -\<V (2Amo)' + « +1 



r+1. s + l 

2 

A-=0 



2 ( - 1 )* [V V + l A* a s A* 13 + M \ r A* a s + * A* 12 

+ 2CV A** « A* 12 { ( MlM2 )^ ( 1 - cos Xl2 ) 



Hftfti) ''*}] P* (COS xia), 



(109) J(ns 1)Xl2 ) = !7r s (2/Am )'- 



• + s + 2 



, + 2 +a [ r+8 A* 13 *A* ]3 +f +1 A* ]2 » +1 A* 12 + r A* ia ' +a A* ls 

& = 



12 



-4^C a B ('■ +1 AVA* 18 +'AV +1 A* 12 ) (1-cos Xl2 ) 
+ ^ 2 2 CVA'\./A k la (1 -cos Xl2 ) 2 ] P, (cos Xl3 ), 



(110) J(r 2 s 1 , Xv2 )=iy(2hm Y 



+ .9+2 



r+ T 2 (-lf [> 12 ;+s A* 81 'A* 12 +r +1 A* sl ' +1 A* ia +^'-A* 81 *+ i, A* ia 

+ 2C 2 K W' w WyM 2 W /2, A.V +1 A* 12 ) {2 ( Mi m 2 ) V2 (1 -cos Xl2 ) 
+ CVAVA\ 2 {2 Um 2 ) Vs (1~cos xi^-UMa)- 1 ''} 8 ] P. (cos Xl2 ). 



(M1M2) 1/i! } 
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The Integration with respect to p. 

§ 7 (F) On referring back to equations (81) to (88) it is clear that before executing 
the integration with respect to p, in the calculation of b (rs) and c (rs), we must 
subtract from I(?-s, X12) and J (rs, X12), as gi yen by (107) to (110), their values 
corresponding to X12 = 0- Now when X12 — 0? we have 

1 -COS X 12 = 0, P A (COS X12) =* 1. 

Thus where P*(cosxi2) occurs alone in (107) to (HO), it must be replaced by 
P A (cos X12) — 1 in the expressions for 6 (rs) and c (rs), the terms P A (cos X12) (l — cos X12) 
and P/, (cos X12) (l-^-cos X12) 2 remaining unchanged, since the corresponding terms in 
I (rs, 0) and J (rs, 0) vanish. 

The variable p is involved in b (rs) and c (rs) only through p dp and xn> the latter 
being also a function of C R . We may therefore formally execute the integration 
with respect to p by writing 

/.CO 

(111) 0*i2 (C„) = (2A-+1) (AHMa)~ Vi C B {l-P*(cOSXi 8 )} J3#, 

Jo 



»oc 



(112) 0'* ]2 (C n ) = (2&+1) (miM 2 )"' /2 C b (1-cosxia) Pt(cosxi 2 )i>^, 







,00 



(113) <i>"\ 2 (C B ) = {2k+ 1) ( MlMa )- 1 /. C B (1 -cos Xl2 ) 2 P* (cos Xl2 ) p #. 

Jo 

The nature of these functions depends on the law of inter-action between molecules 
at collisions, and by keeping this law unspecified we retain the utmost generality 
in our theory, which implies no property of the molecules save that of spherical 
symmetry. 

By means of the well-known equation 

(114) (*+ 1) P*+i (cos X )-(2A+ 1) cos x P ft (cos x ) + kP k _ 1 (cos x ) = 

the function cj> k l2 (C E ) can be expressed in terms of $\ 2 (C K ), for different values of h, 
as follows : — 

(115) ^ 12 (cy = J±l ^\ 2 (c B )-^ 12 (o B ) + ^^-^(Qh), 

and by a repeated application of (114) we may obtain a similar expression (involving 
<j> l l2 (C R ) for I = k, k±l, h±2) for f\ 2 (C B ). 

To avoid unnecessary formulae, we shall not write down the forms taken by 6 (rs) 
and c (rs) on substitution of the results of this section till after we have considered 
the. next step in the integration. 
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The Integration with respect to C and C u . 

§ 7 (G) In the expressions for b (rs) and c (rs), integrated out with respect to all the 
variables save C and K , it is now convenient to make the transformation 



(116) 



x 2 = hm G 2 9 y 2 = Am C u 2 . 



In connection with this we shall use the following notation :• 



(117) B* iaa (m, n) = (2Am )» + ^A* ia »A* 12 = (2Am )* + ».*A*( A i 1 C a , M). n ± k (fiiC 2 , M) 

= *A* (2 Ml /im C 2 , 2^ 2 hm G R 2 ) . n A* (2^hm G 2 9 2vJim<P*)> 
= m A* (2^, 2/x 2 2/ 2 ) . B A* (2^, 2/xgy 2 ), 

(118) B* m2 (m, n) - (2Am ) BI+tt .^A^^A^, 

= m A* (2/x 2 ^ 2 ? 2 Ml? / 2 ) . »A* (2/unx 2 , 2fi 2 y 2 \ 

We have here used the fact — cf. (96)-(l00) — that n A. k (p 2 9 a- 2 ) is a homogeneous 
polynomial of degree 2n in />, o% 

We now use equations (.83), (84), (107), (108), in conjunction with §7 (F), to 
write down the following expressions* for b (r, 5), taking particular note of the signs 
of the various terms :— - 



r + l, s + l 
ft = 



(119) 6 ia (n.s,) = - 1 /- Vl v, e-*+* ± [>* 13 {r 12 y) {B* (r+ 1, s) + B* (r, s + 1)}_ 



+ 4,u.#y*i2 (risy) B* (r, s)] 1212 »y cfe cfy 



(120) & 13 fo*i) = 



1 f* 



e-<*-^ ,+I is +1 (-l)*[0* 13 (T ]2 y.)W / 'B*(r+l, S ) 



7c = 



-2 (^i^-'VB* (r, s) WB* (r, s + l)} 

- 4 (miM 2 )' V<A 2 (ti 2 2/) B*(r, 6')] 2112 xS/dxdy. 



In a similar way, from (87),- 88), (109), (110) we obtain the following expressions 
for c (r, s) : — 



(121) c 12 (r 1 s 1 ) 



1^ 2 | e"^ S [0* ia (T^){B*(r + 2, s)+fB*(r+l, s + l) 

+ B* (r, s + 2)} + 8/x^V^ (t^){B* (r + 1, 5) ' 

+ B /c (r ? 5"+ 1)} -16j* 2 2 y i </>" k i2(Ti2y)B k (r 9 s)] m2 x 2 y 2 dx dy. 



* In (119)-(122) the 'suffixes 1212 or 2112, which should be appended to the symbols B k (m, n) — the 
same for all those within any one square bracket — are for convenience of printing indicated only by 
being placed after the bracket itself. 
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r+2, S+2 



(122) cuirfr) = % VlVa \ e-^ 2) 2 {-l) k \j>\,{r l2 y) {^W(r+2, *) + $B* (**+ 1, s+l) 



k=0 



+ A* 21 B* (r, 8 + 2)-4(/x lM2 )- 1/2 2/ 2 Ui^'B* (r+ 1, s) 
+ M2i 1/9 B*(r, ^+l)j + 4 MlM2 y 4 B^^ s)} 

-2 (^-VB'fa *) W''B*(n s + l)} 

- 16mi//. 2 ?/V //c i2 (tis^) B* (r, s)] 2112 aY dte dy. 



The symbol t in <56\ 2 (t 12 ?/) is defined by the equation 



(123) 



T12 



^ 



m r 



The integration with respect to x in the above expressions is of a quite elementary 
nature, but it will not be executed ki general terms owing to the complexity of the 
polynomials W (r, s), which are integral in x 2 . Any individual term in the integrands 
of (119) to (122) is of the form (so far as concerns x t 



(124) 



-CO 

6 -*V (w+1) dx = iir 1,9 (m + ^) m . 

Jo 



The integration with respect to y will similarly not be executed in general ; in any 
case, owing to the unspecified functions <f> k 12 (ry), this integration could be only 
formally completed, and until we come to consider special types of molecular models 
we shall be content to leave b (r, s) and c (r, s) in the above form. 



The Complete Expression for A 12 Qi. 

§ 7 (H) On referring back to § 7 (A), and the definition of b (r, s) 9 c (r, s), it is clear 
that we are now able to write down the complete expressions for A 12 Qi in the two cases 
we have considered. This involves taking into account all the terms (r = to 00 ) in 
F(U, V, W), with their appropriate coefficients, as in (78), (80); and in order to 
make the expressions more symmetrical, it is convenient to change the values of Q x 
slightly, by multiplying them by certain numerical factors (ef. 26, 30). Thus 
writing 

(125) ^ ~ X 



X 



r— l,s— 1 — 



1.3.5 ...(2r + 3)r.l.3.5...(2* + 3)« 



(126) 



^« = 



1.3.5... (2r + 5) .1.3.5... (2s + 5)' 
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we have the following equations for A 12 Q X in the two cases under consideration : — ■ 

(127) i-Jf—^^^fi^ = B i^2 o X 1 .„ ] ; s _JA.-M^(^^)+M 12 ,/2 A.- ] , 2 ^ 2 (^i)} 



{2hm 1 ) s+1 



QO 

AtflftC! 2 * = C Cn 2 X'„{y r , I c 12 (r- 1 s 1 ) + y r , 2 c 12 (r 3 « I )}. 



l*=0 



^ ' 1 .3.5..*(2s + 5) 

The corresponding values of AnQi in the two cases are obtainable from (127), (128) 
by replacing the suffix 2 by 1 throughout. We will write b n (riSi) and c n (r^i) 
respectively for the values taken by 6 12 0"i$i) + &i2.0Vi) & n d c l2 (r^-t-c^ (r^) when the 
distinction between the suffixes 2 and 1 in these expressions is abolished. In place of 
//i and jul 2 we now write |-, and /x 12 — //^ = 1 , m = 2m 1? while W 12 (m, n) and 
W 2l (m, ?i) become identical, and equal to 

(129) m k k (9i?,tf) n & k {Q(?,i?) = W(m 9 n). 

It is convenient to express &n (n^i) and c xl (r^i) in terms of c/> 7t n (rt/) only, eliminating 
(p fk ii{ry) and //A n (ry) by means of (115) and a similar equation for <j>" k ii(ry). When 
this is done it is found that the coefficient of <j> k u (ry) vanishes for odd values of k, on 
account of the factor ( — l) 7 ' in fr^^Vi) and c^^Vi)- The following are the results 
thus obtained* : — • 



(130) Mn*i)-W 



0\ s) 



«/ *j 



k= 1 



e _b. +y ., 2 ^ ^ Tu ^ B% ( r + j ^ g ^ + B24 ^ s + x ^ 



2y 



J2k+1 



2k 



U& + 



-B 2 * +1 (r, s) + -p- B 2 "- 1 (r, s)-B 2k (r, s) 



as ^/ cfe c%. 



(131) d, (r lSl ) = J ^ 2 e- fe3 ^ 2) 2 f k u (my) 



[»\ s] 



7; = 1 



B 2 * (r + 2, «) + |B 2/< (r + 1 , s + 1 ) + B 2 * (r, s + 2 ) 



+ 4y 2 |||±i (b 2 " + x (r + 1 , «) + B 2 * + * (r, s+l)\ 



+ ■ 



2& 



4^; + 



_/B 2 *- 1 (r+l,s) + B 2 *- 1 (r,s + l))-(B 2/c (r+l,s) + B 2/i (r,s+l)) 



+ 



j(2/c + 2)(2^ + l) B2 , +2( } 
J 1 Uk+3 )Uh+\) v ; 



(4& + 3)(4fc + l) 
(2^+l) 2 



+ 



+ 



44 s 



l (4^ + 3)(4/c+l) (4&+l)(4£-l) 

2k(2k-l) m _ 2( , 
(ik + l)(4k-l) { ' 

2 (|iB«(,,) + Ji I B-(, s) 



+ i B 2 * (r, s) 



xy dx dy. 



* In these expressions <£*n (rny) is the equivalent, for an encounter between two molecules of the same 
kind, of <£ & i2(Tny) for molecules of different kinds. Thus (cf. Ill) 



rao 



^n (niy) = 2 (2& + 1) R {1 - P* (cos X n)}i^ 



o 



where the law connecting xn with p and C^ may differ from that for xi2- Also r n now becomes (2/im 3 ) Va. 
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The term corresponding to k = is absent in both the above cases, since 0° n (my) 
is itself zero, so that (j> 2 u (my) is the function <p n (my) of lowest order (k = 2) in 
&n (n^i) or en (nsi). The upper limit of & in the case of &n (n^i) is equal to the 
integral part of the lesser of the two quantities J (r + 1) and %(s + l); this is denoted 
by (r, ,§). Similarly the upper limit of & in the case of c u (nsi) is the integral part 
of the lesser of the two quantities |-(r + 2) and ■% ($ + %), which we denote by [r, s]. 
Thus, when r = or s = 0, &n (n^i) = 0. 

We can now write down the complete expressions for AQ, in the two cases above, 
as follows : — 

( 132 ) (f™ 1 ?* 1 oX AUxC^ 
1.3.5 ...(2s + 3)s * x 

1 r^T °° 

= B i |i 2 X,_m-i [&-i, 1 {&11 (n»i) + b 12 (nsx)} +m 1 2 1/2 /3 1 .- 1 , 2 &i 2 (r^i)], 

( 133 ) 1 q^^IT I r\ Au ^ 25 = C ° c " 2 X , „[y,.i{cn(nsi)+c ls (r 1 s,)} + y r . s c ia (r A )]. 
1.3,5... (2s + 5) c=o 

In the present paper we are concerned with the application of these formulae only 
to simple gases, in which v 2 = and hence b }2 (r^j) = b 12 (r^Si) = c 12 (r^) = Cj 2 (?Vi) = 0. 
It is convenient to write the reduced equations in the following form : — 

(134) {2hm) s+ l § A UC 2{s+1) = ^ P 2 && w 

(135) (2Mp^ 45 A 2Cps = |, 

^ 7 1.3. 5. ..(2«+5) 2* - 



r ~ 



In (134) we have substituted r+l, s+1 for r, s in (132), multiplied by 3/v, and 
used the notation given by 

(136) b rt = -B^ t9 b n (r + l 9 s + l); 

the first term in (132), with factor /3_i, vanishes, since &n(0, s) = 0. Similarly in 
(135) we have written 

(137) * c rs = — CoX^Oii (n*0- 

§ 8. The Expressions for the Coefficients in the Velocity- 
Distribution Functions. 

§ 8 (A) We have now obtained expressions for AQ, the rate of change of a function 
of the molecular velocities due to encounters, in two different ways : in § 3 AQ was 
found from the equation of transfer, while in §§ 4-7 it has been determined by direct 
calculation, By comparison of (26) and (134) — substituting 5+1 for s in the former 
— and of (30), (72) and (135), we deduce from these different expressions for AQ the 

vol. ccxvi. — A. ^ 2 u 
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following linear equations connecting the unknown coefficients in the velocity- 
distribution function : — 

(138) 



2 (3 r b rs = 1, 

r = 



(139) 



oo 



*Yt TS * 

r = 



These are true for all values of s from to oo } the coefficients b rs and c rs being 
completely determined, in terms of the molecular data, by (130), (l3l), # and (136), 
(137). If we assume that certain convergency conditions are satisfied (138) and (139) 
lead (in the way usual in the case of a finite system of linear equations) to the following 
expressions for f} r and y r : — 



(140) 



ft = V, (6„)/V (&„), y r = V r (C rs )/V (c n ), 



where V (b rs ) and V (c rs ) denote the infinite determinants formed from the arrays (&„) 
and (c„), thus, 



(141) 


v (b rs ) = 


K 


&10 


"20 


^30 






K 


&U 


&» 


31 






K 


&18 


&22 


t>32 • 






&08 


&13 


&23 


^33 



c oo 


^10 


&20 


C 30 


C 01 


Cll 


C 21 


c 31 


C 02 


C 12 


G 22 


C 32 


C 03 


c ys 


°23 


c 33 



and V r (6 ra ), V r (c„) denote, the determinants obtained by replacing each element of 
column (r) in V (b rs ) and V (c rs ) respectively by unity. 



The General Expression for the Velocity Distribution Function. 

§ 8 (B) This completes our solution of the fundamental problem of this paper, i.e., 
the determination of the velocity-distribution function for a " nearly perfect " simple 
gas, composed of monatomic molecules of the most general type, and which is slightly 
non-uniform as regards temperature and mass-velocity. The solution will be sum- 
marized as follows (cf. (10), (73)) : — 



(142) 



f(\J V W) = ( m ) *0-^ m (u 2 +v' 2 +w' 2 ) 

7T / 



(2 Am) 1 



l ~~ B °Y\ U te + V Yy + W TzJ^o US.5...(2r + S)r 

C (2hm) (c 11 U 2 + c, 2 V 2 + c 3 3W 2 + c 2 3VW + c 31 WU + c 12 UV) 



A._iG 



\2r 



oo 

2 



(2Am) r 



= o 1.3.5...(2r+5) yr 



C 



2r 



The suffix 1 throughout these equations may now be omitted. 
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where c n , c 12 , &c, are given by (72). The coefficients /3 r and y„ for r = to r = oo, 
are given by (140), where (cf. (136), (137), (130), (131)), 



(143) b rs 



32B „A„ \\e~ v ^ sV*(Ty) 

k = l 



B 2/c (r+2,s+l) + B 2 *(r + l,s + 2) 



+ 2?/ 2 ||l ± ^ B 3 ** 1 (r + 1,8+1)+ -^- B*- 1 ^ l,a+l) 
[4&+1 ' 4&+1 x ' 



B 2k {r + l,s + l) 



oc y ctrju cox/ } 



(144) c„ = 72G v\' 



[r, s] 



e-^'W) s <t>*{ry) B 2 *(r + 2,s) + fB 2 *(r+l,s + l) + B 2 *(r,s + 2) 

ft = l 



2 | 2A + 1 



+ 4?/ 2 ||g^(B a+1 (r + l,s) + B a+1 (ns+l) 



-u/v 



+ 4f ± Y( BM "V+l > *) + B 2i -V } s-l))- [B 2k {r + l,s) + B 2k (r,s+l) 



f f(2* + 2)(2*+lW ( } 
^ \(4&+3)(4£ + l) v ; 



+ 



(2&+l) ! 






4# 



(4& + 3) (4& + 1) (4* + l)(4A-l) 

, 2/c (2a;— 1) r)2s-2 /„ „\ 
+ (4&+l)(4&-l) K ' } 

_ 2 ( 2 4±1 B»+i ( r . s) + -4^-r B 2 *- 1 (r, s) 



+ l)B 2 *(r,s) 



Ah+ 



1° K,S)+ ih+l 



xy dx dy 



where (cf. (123)) 

(145) 

and, by (ill), 
(146) 



= (2hm)' 



.00 



<p 2 " (z) = 2(U+l)z\ { 1 -P a (cos x )}p dp, 



where P* (cos x) is the usual Legendre's coefficient, and x ^ a function of p and z 
which depends (§ 4 (D)) on the law of inter-action between two molecules at an 
encounter. The factors \ rs and \'„ are defined by (125), (126), while the functions 
B* (n s), which are integral polynomials in x and y, with merely numerical coefficients, 
are defined by (129) and (96). In the upper limit of jfe, [r 9 s] denotes the integral 
part of the lesser of the two quantities j%r + 1 , ^s + 1 . 

The factors B , C are, as yet, arbitrary; we now assign to them the values 
determined by the equations 



(147) 



'00 



1 



"00 



This makes B and C each equal to v~ x multiplied into a function of (2hm), i.e., of 

2 u 2 
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the absolute temperature. The elements &„, c rs and the coefficients /3 r , y r then become 
fun ctions of the temperature only. 

The coefficient j3^ is not determined by the above equations, but is given (cf (75)) by 



CO 



(148) 5_!= - 2fr/(r+l). 



f = 



Properties of the Determinants V (b rs ), V (c rs ). 

§ 8 (C) On inspection of (143) and (144) it is evident that 

^I4yj o rs = o sn e rs c 5f , 

so that V (h rs ) and V (c rs ) are symmetrical determinants. 

In expression (143) for b rs , the variables of integration, x and y, are never negative, 
so that (cf. (129), (96)) r A*, *A k and B* (r, s) are essentially positive (or zero) for all 
integral values of r, s, and h ; further, since P /c (cos x) never exceeds unity, <j> 2k (ry) is 
also always positive. It is evident, therefore, that b rs must be essentially positive 
if this can be proved true of 

(150) B 2 *(r+2, s + l)-2y 2 B 2k (r+l, s + l) + B 2 *(r + l, s + 2). 

Now 

(151) B 2 *(r + 2, s + l)-y 2 B 2k (r+l, s+l) = s+l A 2k [ r+2 A 2k -y 2r+1 A 2k ] 

.« = at (*+£)* (*-2*)t ^ t = sk(t+$) t (t + 2k)\ ■ 

^.s+l\2k(y\ ' '4 (r+ l)^_i {r + ^r) t _2k-l „2t 7 .2(r+2^t) f /„ , Q ;\ /; Q/.N , f / ■ 3\] 

~ W «r a (*+i), (*-2A)i xy {{r+J t){t u > +t ( r+i *»> 

every term of which is positive. Interchanging r and 5 in (151), and adding the result 
to the latter, we obtain (150), which, with b rs also, in consequence, is essentially 
positive. 

From (151), moreover, it is clear that the numerical coefficients in (151) or (150) 
increase with r or s, and the same is readily seen to hold good also in the case of 
B k (r> s). As r or s increases, therefore, the numerical coefficients and the degree 
(in x and y) of the integrand of (143) increase, while if both r and s increase, new 
positive terms are added to the integrand. Hence, provided that the functions 
(p 2k (ry) satisfy certain simple conditions, # b rs steadily increases with r or s, and the 
consideration of even a single term of (151) or the integrand of (143) shows that this 
increase is without limit, Le.> b rs tends steadily to infinity with r or s. 

* It is easy to see that the increase with y of <f> 2k (ry) is less rapid than that of y ; if <j> 2k (ry) is constant 
or steadily increases, though less rapidly than y, b rs will steadily increase with r or s. But much less 
restrictive conditions might be devised, e.g., if <j> 2Jc (ry) decreases like y~\ the above result still holds good. 



$+l\2kty 
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I have little doubt that, with rather more trouble, c rs could be shown to share the 
above properties of b rs , but I have not made any serious attempt to prove this ; from 
the numerical calculations in § 10 (A) it appears probable that the increase of c rs with 
r, « is more rapid than that of b rs . 

Properties of the First Row or Column of V (b rs ) and V (c„). 

§ 8 (D) The numerical values of b rs and c rs obtained in § 10 suggest that many 
further general properties of these elements might be determined, with sufficient 
trouble, and that the convergence of the determinants V (b rs ) and V (e rs ) might thus 
be demonstrated. Owing to the considerable algebraic difficulties involved, however, 
I have so far made little progress towards the proof of such properties, except for 
the case when r or s is zero, i.e., for the elements of the first row or column of 
V (b rs ) and V (c r8 ). It will be shown that 

for all values of r. 

This will be proved as a particular case of the more general result that 

(152 a) (s + 1) b rs (k) = c rs (k) when the lesser of r and s is even, and k = [r, s], 

where b rs (k), e rs (k) denote the parts of b rs and c rs respectively which are due to a 
particular value of k in (143), (144), while [r, s] denotes the upper limit of k, as usual, 
i.e., k = fr+1 or -Js + 1, whichever is the less. Thus if we suppose that r §= s, and 
that s is even, (152a) takes the form 

(153) (s+l)b„{is + l) = Cnds + l). 

When 5 = 0, this value of k is unity, and 5 r0 (l), c r0 (l), which usually form only 
a part of 5 r5 , c rs9 become the whole, so that (152) is the particular case of (153) 
corresponding to this value of s. 

Since B /c (r, s) is zero when either r or s is less than &, some of the terms in 
6«(|"S + l)j 0«(isr 5 + l) vanish. In fact, as may readily be seen from (143), (144), we 
have 

(154) b„ (i* + 1) = B2B v\ rs e- (x "- + '^f +2 ( T y) \ B s+2 (r+1, s + 2) 



*> *> 



+ 2_(i±2) yS B . + i ( r+ i s + i)\ x y dx dy, 
2s + 5 J 



(155) C n (js+l) = 72G v\' rs j| e-^f+ 2 (ry) JB s+2 (r, s + 2) + ^£±|1 ^B' +1 (r, s + l) 

+ ii f + !)(« + 2) y i B s u A a-y ^ ^ 
(2* + 3)(2* + 5)* v ' ; j y y 
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From (96) it is easy to see that 



(156) 
so that 



k\k 



k\ 



(k+i)k 



TV ( x v) k > 



(157) B s+2 (r + l, s + 2) + ?i^±|l 2/ 2 B s+1 (r + l, s + l) 

= ( S + 2 ) ! ^ ( r+1 )t (r + $)t-.-l x 2t 2 (r + .+3-«) 

(« + f) s+2 t = s + l (t + i) 4 (*-S-l)I 

By putting r in place of (r + l) in (157), and adding {2(s + 2)y 2 /(2s + 5)} times a 
similar expression in which r, s replace r+l, s + l in (157), we also have 



(158) 



V ' ; 2s + 5 y V ' 7 (2s + 3)(2s + 5) y v ' ; 



(s + 2)\ ^ 

(s+f) s+2 < = s (t+$) t («-«)! 



r 

2 



r * . V " ^ 27£-s aJ 2^2(r+* + 2-t) 



We now substitute the expressions on the left of (157) and (158) into (154) and 
(155), and integrate with respect to x by means of the well-known formula 



(159) 

we thus obtain the equations 



(160) & f ;(|s + l) = 8B i/A„ e-^' +a (ry) 



,00 

e -*V«+i>da. = iV /2 («+J) t ; 





(s + 2)\ '+. 1 (r+ l) t (r+f )«_,.! 2(f+s+4 _ t) , 
(s+f) s+2 t = .+i (t-s-l)I 



(161) c„(i* + l) = 18(W«X' 



rs 



„-y*.s+2( ,.\ ($ + 2)1 ^ 9j V +yX~s ,,2(r+s+3-~a) ,7 ?/ 

6 * M (^|L 2 i («-,)! y ^ 



or, changing the notation so as to make the lower limit of t zero, and inserting the 
values of X m \' rs according to (125), (126), i.e., 



(162) 



we have 



X 






(r+s+i) 



rs 



(r+l)(s + l)(r+f),. +2 (s + f) s+2 ' 



A' 



2 



■ (r+ff+4) 



(^ + l)r + 2(s + l) S + 2 



(163) 6„(i*+l) 

= 2- (r+S+1) Bn^7T I/a 



(164) c„(i* + l) 



r. 



(* + 2)l 



( S + l)(r+|l +2 {( S +f) s+2 } 2 



r— s 



e~^' +a M s r-A(r+f)«^ (,+3 - <) ^, 



< = 



<y» 



9 . 2- (r+s+3) CW 3 ^— i? 



(8 + 2)1 



(^ + f)r + 2 {(s + f) S+2 } 2 J 



r—s 



e-^'(ry) 2 ^(r + f)^ 3 -"^. 



t = 



THEORY OF VISCOSITY AND THERMAL CONDUCTION, IN A MONATOMIC GAS. 3 1 7 
The ratio of these two expressions is given by 

rs Xl&S i ->- ) -"L 

Crskks+l) ~~ S+l 9 

where A is a quantity independent of r and s. When s = 0, as we have seen, 
b rs (t^s + 1 ) and c rs (j^s + 1 ) become identical with b rs and c rs respectively. Hence 



h 

u r0 



'r0 



xx, 



and since B and C have been chosen so that b 00 = 1, c 00 = 1, the value of A must be 
unity. Hence, when s is even and r ~ s, 

(5 + 1) b rs (£*+ 1) = c rs (|s + 1), 

with the consequence that 

o r0 = c r0 
as a special case. 

It is convenient to introduce the notation 

fee 

(165) e-»V» M V 2im+1) dy = 1^ (m + li^-sw, 

Jo 

so that if ^{ry) had the value unity, the value of K m _ 2 *,* would also be unity, by 
(159). In terms of this notation (163) and (164) may be written as follows : — 

(167) c r5 (**+ 1) = 9 . 2-<— 5 > Co^rr, r/l^f \ a S ,„AK r _ M)ts+1 , 

By writing £ = r—s— £' it is evident that 



r-s r~s 



£=0 «=o 

By giving to r and s in (166), (167) the value zero, we have 

(169) & 00 = -g-f 5 B J/7rK 0) u C 00 = ^CoWrKou 

whence, remembering that (c/1 (147)) B and C are so defined as to make b m and c { 
each equal to unity, we have 

(170) -t>o — ~4^ t? , ^0 — 25 ^ , B = ^O . 
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We will substitute these values into (166), (167), and write s = 0, in order to obtain 
expressions for b r0 and c r0 as follows : — 

where we have written 

Kkf/i) K r "ZZZ £i Zj r \j£ jS-i J, 

* = 

so that 



"oo 



(173) *, = K , 1 = -^ j; e-»y( T y)y«dy. 



19, 







It is of interest to examine also the r th successive difference of b r0 or c r0 , which we 
shall denote by S r0 b rQ or § r0 c r0 . We have 

(174) § r J(r) =/( r )- r C 1 /(r-l) + r C 2 /(»-2)- ... • 

Then, from (172), it is easy to see that 



r-m 



d r QK r — A Zj \ Lj r V7 m Zj r „ m V7j { J\. ij 



m-0 t-0 

r r—t 

(-9,) m .f! 

711 



— 2 ' 2j r Kj t &. tl Zj ( — 2) m r _ t Kj 

t = Q m = 



— ( — 2) r 2 ( — l)* r C t K iil3 



« = 



since 



Hence 



i £ (-2)v.c w = (i-2)'-« = (-i) r -*. 

m = 



(175) <5Ao = i^ = (-2)- r (^o)- 1 2 (-iy r G t K tA . 

t = Q 

Similarly the r th difference of 



r—s 



2 r s 2u r _ s Kj t K. ti ii 2S+u 



which is the part of b rs (^s+l) or c n (^s + l) which depends on r (s being even and 
r == s) is equal to 



rv* ^ & 



vV + s ™ 2 ( — 1) ) ._ s O i K i) i/ 2Hl( 



00 CO 



Symmetrical Expressions for 2 f$ r and 2 y r . 



r = r = 



§ 8 (E) While V (b rs ) and V (c rs ) are symmetrical, the derived determinants 
V r (b r8 ), V r (c„) are necessarily lacking in symmetry, and our expressions for )3 r and y r , 
when we attempt to make successive numerical approximations to their values 



THEOEY OF VISCOSITY AND THERMAL CONDUCTION, IN A MON ATOMIC GAS. 319 

for particular types of molecular models (cf. § 10), appear not to converge at all 
rapidly. Fortunately, in our applications of the velocity-distribution function to the 
theory of viscosity and thermal conductivity, we need to know not the individual 

00 00 

values of the /3's and y s, but only the sums 2 /3 r and 2 y r ; for these it is possible 

r — r = 

to determine .symmetrical expressions which are found, in practice, to be highly 
convergent. 

In what follows we shall use the symbol § mn placed before a function of the 
integral variables r 9 s (such as b rs or v rs ) to denote the (m, n) th successive difference 
of this function with respect to r and s respectively. Thus 

<W(n s) =/(r, ^-^/(r-l, s) + w C 3 /(r-2, s)- ... 
$o*f(r,s) = f(r 9 s)- n Cif(r,s-i)+ n C 2 f(r 9 s-2)-... 

When we substitute 6^ or c„ for /(r, 5) in the above formal expressions, any term 
with a negative suffix is to be omitted as being zero. 

Since the value of a determinant is unaltered by subtracting from the elements of 
any one row or column the corresponding elements of any other row or column, and 
since this process can be repeated indefinitely often, it is clear that from ( 141), by 
subtracting the ($ — l) th row from the s th , for all values of s from 1 onwards, we have # 

(176) V (&„) = V (S 01 b rs ), V (c„) = V (S ol c„). 

The same process applied to V r (6 rs ), V r (e rs ) leads to determinants identical with 

V (§ 0l h rs ) and V (S Q1 c r8 ) respectively, save that in the r th column -all the elements are 
zero except the one in the first row, which is unity. Evidently, therefore, V r (b rs ) 
and V r (c rs ) are the r th minors of determinants which are respectively identical with 

V (S Q1 b r8 ) and V (S 0l c rs ), except that in each case all the elements of the first row have 

CO CO * 

the value unity. Consequently the sums 2 V r (b r8 ) and 2 V r (c rs ) are equal to the 

sums of the minors of the two determinants just described, i.e., they are equal to 
these determinants themselves. Thus, by (140), 

{•\>7*7\ V /Q — Vgoig rs) -y ._ Y_ \VQl rs) 

11771 r-0^~V(^6,.)* -0 7r V(V«)' 

where we have 

(178) &' ,= 1, c / ,= l, (s=Otooo) ? &'„ = &„, c'„ = c r „ (^ = 0tooo ? r = ltoco). 

* When s = 0, <S 01 should be replaced by <$ 00 . 

vol. coxvi. — a. 2 x 
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To the above determinants we now apply the same process of differencing by columns 
which has already been applied by rows, and we thus obtain the equations 

K ' rto Pr ~ V(5 n 6„) ' rt 7r ~ V (^ llC „) ' 

The determinants V (S u b f rs ) and V (S n b rs ) are identical save in their first rows ; all the 
elements of the first row of the former are zero save the first, which is unity. Hence 
V($ub' n ) is equal to the principal minor of V (S u b rs ) ; we shall denote it by 
^' ($iJ>rs)' Hence, and with a similar notation for the principal minor of V (§ n c rs ), 
we have # 

I B = v'M„) 2 v - YLfe.) 

,t ^ V ( J u &„) ' rt ^ ~ V (j uC J 

All these determinants have now regained a symmetrical form. 

It is convenient, partly for the sake of elegance, and also because it imparts a 
highly convergent form to the elements of our determinants (c/! § 10) to continue 
this process of differencing still further, as follows. We repeat the whole of the 
above operation of differencing by rows and columns an indefinite number of times, 
beginning now at the second row and column (thus leaving unchanged the values 
both of V and its principal minor), and afterwards successively at the next later row 
and column than on the previous occasion. The general element thus becomes S rs b rs 
or .<Lo, c , and we have 



'rs^rs") 



(l80) 2 8 ~ K^rs^rs) 2 V = \°rs C rs) 

r=0 T V (§ rs b rs ) ? r==0 V ^ (S rs C rs ) ' 

where the dash (') denotes the principal minor of the corresponding determinant. 
These expressions could, of course, have been obtained directly by a re- arrangement 
of the original equations of transfer, but it seems preferable to use the latter in the 
more simple, natural forms chosen, and to make this transformation by differencing 
in relation to the determinants formed by the elements b rs , c rs . 



§ 9. Consideration of Particular Molecular Models. 

§ 9 (A) While, as we have seen, certain general properties of the elements b rsy c rs can 
be demonstrated without the assumption of any property of the molecules save 
spherical symmetry, it is possible to carry our investigations much further when we 
represent the molecules by particular models of simple type, such as point centres of 
force, or rigid elastic spheres. This involves, primarily, the examination of the 
functions cj> 2k (ry). 

* When r or s is zero, the corresponding suffix of S n should also be written as zero. 
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Molecules which are Point Centres of Force varying as r~~ n . 

§ 9 (B) When the molecules are point centres of force varying inversely as the n th 
power of the distance, the angle x in the expression (146) — -cf § 4 (D) — is given by 
the following integraP : — 

(181) x = 2 [W-O,/*)"- 1 ]- 1 ''^. 

Jo 

Here % is the least positive root of the equation [1 — if— (^/a)"" 1 ] = 0, and a is a 
multiple of p, thus, 

where K is a constant which measures the intensity of force between two molecules 
at unit distance. Hence (cf (ill)) 

2 

(183) 4>*{Ty) = 2(a+l)(^-Y 1 C u 1 -^-i C {l-P 2fc (cosx)}ad* 

\n-—l/ Jo 



2 ^^)fer(£f 5 i>- p "(-x)}^ 



A ( y 2 Y 01 "" 1 * 



\2hmJ 

where n A k is a constant depending on n and k, but not on y or h (i.e., not on the 
absolute temperature). 

When this value of <p 2k (ry) is substituted in our expressions for b rs and* c rs , it 
becomes possible to execute the integration with respect both to x and to y in terms 
of gamma-functions. Thus (cf (165)) 



.00 



(184) e-»y* M f ( " + » dy = JV'' (m + $) m K w _ 2 ,, k 



n^-^-k 



_ * / 1 \ 2{n ~ J) 

2hmJ 

n~5 



oo / , n— 5 



/ , n-b \ 

_,.2 2(m+l+ 5 — :) 7 

e v y \ n lJ ay 







iA biz: r(m+2- 



K 2hmJ \ n — 1/ 

so that • 

n — 5 t 

(185) , = Ko,, = 1 b 6T -v. ii a 1 ^)' ( "" 1) r(4-^ T 

* Qf. § 14, p. 454, of my former memoir, * Phil. Trans.,' A, vol. 211 (1911). The V of the formula 
there given is the relative velocity of two molecules, which in our notation is (/xi/x 2 )~ 1/2 C R .== 2C R when the 
gas is simple. 

2-v O 
JSl. Zi 
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to ~ 5 , x U + 3 — 



(186) K tl = 2 ^,A,(^-y-» {((+f ) M }-r( t+ 4- ; A I ) = V ((+ ^ 1 " K>,. 

Hence the values assumed by B and C in this special case are as follows {of. 
(170), § 8(D)):- 

(187) B = 153 (2Am) 1/2 ^T c ^375 (2hmf^ 

32 7 r\A i r(4-- 2 — V ° 8 T\A,r(4 2 



/ \ n — 1. 

From (186) we have 

r 
O rQ = C r0 = /C r //C = ^j 2; y^Ji.^ 1 /iV 0t J 



i = 

2 



£ + 3- 



= 2~ r 2 ^ V n-lA 

= 2"'F(-r,4 ^-,£, -1 

in the notation of the hypergeometric function. It may hence be shown, without 
much difficulty, that (if n > 5) & r0 and c r0 steadily increase to infinity with r, the 

rate of increase being comparable with that of r n - 1 . 

Since the functions <j> 2k (ry) all depend on h in the same way, it is clear that, with 
the above values of B and C , the elements b rs and c rs and consequently, also, the 
coefficients jS r and y r in the velocity-distribution function for molecules of this type, 
are independent of h, i.e., they are independent of temperature. They are, indeed, 
pure numbers, depending only on the molecular mass and on the force constant of the 
molecules. 

It is of interest to determine the value of the elements § r0 b r0 (or § r0 c r0 ) of the outer 
row or column of V (8 rs b rs ), in this special case. We have, by (175), 



r 



Ko^rO — °rQ C rQ -— ( — 2) 2 ( — 1) r \J t K. t% 1/-&0, 1> 

t = 

t + 3 2 ' 



= (-2)-Z(-l)' ^ . »-'k 



t 

\ 



= (_2)-'F(-r,4 ~,hl), 

\ n—1 j 

in the notation of hypergeometric functions, or, in terms of gamma-functions, 

r(|)r( / r-i+- 2 - r N ) (V-1-+ 



(188) <y> r0 = s r0 c ra = (-2)-' v n \ = (-2)- r V 7— gv Vr 

r H-*) r k=i-* 1 ( } " 
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As r tends to infinity, the last expression tends to zero more quickly that 

and also it is alternately of positive and negative sign, after the first two terms 
(both S 00 b r0 and S 10 b 10 being positive). 

Maxwellian Molecules : n = 5. 

§ 9 (C) It is now easy to see what are the special properties of the fifth-power law 
(n = 5), the law obeyed by the molecules which we term Maxwellian, which enabled 
Maxwell to work out the theory on this hypothesis with such great simplicity and 
accuracy. When n = 5, we have from (186) 

(189) K« il = 6 A 1 , 

which is independent of t. Hence, by (170), (171), (172),' 

(190) K, x = Ko ! = *Ai, k t = 2~ r 2 r C,K f , , = 2~ r K 0jl 2 r G t = K , x = * , ho = <Vo = 1 ; 

t = o t = 

(191) 8 r , b r , = S r , c„ = (-2)-' 2 (-l)',C« = (-2)-'(l-l)' = 0, (r>0); 

t = 

(192) B = *|A (v b A x y\ C = 25 (.A)" 1 - 

From (191) and the equation b 00 = c 00 = 1 we deduce that in this case the principal 
minors of V (§ rs b rs ) and V (S rs c rs ) are equal to these determinants themselves, i.e., 



GO GO 



(193) 2 /J r = 1, 2y f =l (Maxwellian molecules), 



r = r = 



while from (190) it appears that all the elements of the first row and column of 
V (b r8 ) and V (c rs ) are unity. Hence in V r (b rs ) and V r (c r8 ) the first column and 
column (r) are identical, so that we have 

(194) V r (b rs ) = 9 V r (c rs ) = 0, r>0l 

> (Maxwellian molecules), 

(195) V (6 rt ) - V (&„), V (c n ) = V (c rs ) J 
whence also, by (140), we have 

(196) /3 = 1, y =1, ftr—Jr— 0, (r > 0) (Maxwellian molecules), 

and also, by (148), 

(197) fi_i = — A) = — 1 (Maxwellian molecules). 
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In the case of molecules which obey the fifth-power law, therefore, the velocity 
distribution function has the simple finite form (cf. (142)) 

(198) /(U, V, W) = feY'V^^+vW) 

(Maxwellian molecules), 
where C 2 = U 2 +V 2 + W 2 , c u , c 12 , &c, are given by (72), and (cf (183)) 

(*Oo poo 

(199) 8 A, = 1 (iKm) 1 '* { 1 - P 2 (cos x ) } a da = V" ( Km) 1 ' 2 sin 2 x . a da. . 







•00 



Maxwell # calculated the value of the integral j sin 2 x-^^ ? "the forces being 
repulsive, by numerical quadrature, and found that 



•00 

•2 



ir sin' x . a da = 1*3682, 





so that, for repulsive forces proportional to the inverse fifth power of the distance 
(200) 5 A 1 = ^ 1-3682 (Km) 1/a , 

Lit 

where Km 3 is the force between two molecules at unit distance. 



Molecules which are Rigid Elastic Spheres. 

§ 9 (D) We next consider molecules which behave at encounter like rigid elastic 
spheres of radius <r. This particular molecular model has been more used than any 
other, in researches on the kinetic theory, on account of its simplicity and concreteness, 
which aid the imagination in following or constructing " descriptive " theories of 
gaseous phenomena. As regards the analytical development of the theory, also, it is 
probably the simplest case after that of Maxwellian molecules. The difference 
between the two models in this respect is, however, enormous, the rigid elastic 
spherical molecule requiring the infinity of terms /3 r , y r in the velocity-distribution 
function, just as in the case of the most general molecular model. The comparative 
simplicity of the present model lies in the moderately tractable expressions for b rsy c 
to which it leads. Apart from the methods of the present and my former paper, 



rs 



* Maxwell, ' Collected Papers,' ii, p. 42. His constant A 2 equals it sin 2 x . a da in our notation. 

.0 



GO 
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however, it has not been found possible in the past to obtain any close numerical 
accuracy in calculations based on this molecular model, the errors resulting in previous 
theories (although these have been carefully constructed and closely scrutinized) 
ranging from 10 to 50 per cent. (cf § 11 (F))< 
It is readily seen # that in the present case 

(201) x =0.(p>2<t) $mi x =p/2<r (p^2a) 

so that 

pdp — 2<j 2 sin Jx cos ^x ^X ^ °" 2 ^ n X ^X ~ "~°" 2 ^ cos X- 

As_p ranges from to 2<r, x ranges from to 2tv, and — cos x from — 1 to 1. 
Hence (cf (ill)) 



(202) <f> 2k (ry) = 2(4^+l)C R o- 2 {l-P» (cos x )} d cos x 



-I 



= 4(4^+l)(7 2 (2Am)~%, 

since P 2/£ (/m) dju — 0. Hence </> 2k (ry) depends on h only as regards the numerical 

factor (4&+1), and the present case is, analytically, the same as that considered in 
§ 8 (B) if we write (cf (183)) 

(203) ^— J = 1, or n = oo, and n A* = 4 (4A + 1) a- 2 . 

We may therefore quote from the formulse of § 8 (B) as follows without further 
discussion ; — 

(204) 1V /2 (m + |) m K m „ 2 „ = 2 (4* + l) rr 2 (2Am)~ 1/2 r (m + 2) 

= 2(4^+l)(m + l)!(T 2 (2^m)- 1/2 , 

(205) * = Ko §1 = 647T- V (2Am)" 1/f > K m = K ci (* + 3)«/(*+-|),. 

/onfrt B - 225 (jM! 2 n _ 25 (2M 1/2 



, 2 * 0—64 i/ 2 2 



(207) 6,o = c r0 = 2- 2 ff J4J# = 2-T (-r, 4, f , -l), 

(208) ^Ao = *«fi« = (- 2 y r jy^f r =-(-^\ 2 r-.l)(2r+l)(2r+S)(2r+5)- 

As in the case of molecules which are point centres of force varying inversely as 

* Cf §13, p. 453, of my former paper, 'Phil. Trans./ A, 211 (1911). 
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the n th power of the distance, the elements b rs and c rs > and the coefficients /3 r and y n 
are pure numbers, independent of the temperature (i.e., of h). 

Molecules which are Rigid Elastic Spheres which Exert Attractive Forces. 

§ 9 (E) Experiments on the phenomena of actual gases, as, for example, on the 
variation of viscosity with temperature, indicate that none of the molecular models 
so far discussed in this chapter gives a really adequate representation of these 
phenomena. The best of all the simple models which have been used in the kinetic 
theory seems to be that considered by van deb, Waals and Sutherland, viz., a 
rigid elastic sphere surrounded by a weak field of attractive force. This agrees with 
the known fact of slight cohesion in gases. The effect of this field of force on the velocity- 
distribution function, or on viscosity and thermal conductivity, may be referred 
mainly to the deflections of molecular paths for which it is responsible indirectly, 
through the collisions which it induces between molecules which would otherwise 
pass one another without mutual inter-action, rather than to its direct effect in the 
absence of collisions. The latter effect w T ill be expressly neglected in our calculations, 
which will therefore be inapplicable to vapours in which the cohesion is large enough 
to render this neglect invalid. 

A detailed account of the dynamics of collisions in these circumstances is given in 
§ 1 5 of my former paper, from which the following results are quoted. If the potential 
of the force between two molecules in contact be denoted by f& 2 (reckoning this 
potential as zero when the separation is infinite), the condition that a collision may 
take place is 

(209) P <Po where p = 2a (1 -f g( p ) 

(since the relative velocity, in my former paper written V , is- here denoted by 2G K ). 
The angle x corresponding to such a collision is given by 

(210) . sin J x = p/Po- 

The angle x corresponding to larger values of p, which do not correspond to actual 
collisions, is given by (18 1) if the molecular forces obey the n th power law, but we 
will here make no assumption on this point, as the deflections produced by the inter- 
molecular forces alone will be rejected after equation (211). Consequently 



.00 



(211) f k (ry) = 2(ik+l)cr 2 (l+b 2 /4C n 2 )G n + 2(U+l)G ll \ {l-P„ (cos x )} p dp 

= 4 (4k + 1 ) a- 2 (2hm)S ( 1 + 2h<mb 2 /4y 2 ) +f 2!c (y) 
by analogy with (202) and (183). The latter term f 2/l (y) represents the negligible 



THEORY OF VISCOSITY AND THERMAL CONDUCTION, IN A MONATOMIC GAS. 327 

deflecting effect of the forces alone, as above mentioned, and it will be omitted 
henceforward. Hence, corresponding to the equation (204) of § 9 (D), we have 

(212) K''(< + i). £,_„,* = 2 (4& + 1) , ! (2k)-''T(* + 2) (l+ 2hmb 



4(* + l). 
3S 



= 2 (4k + l)(t + l)lo- 2 (2hm)-' h (l+ j . , , 



where we have written 

(213) S-£|. 

Similarly we have 

(214) k = K 0il = 64ir- V (2hm)- lk (l + '-§,, 

/o lfl \ -R _ 225 ( 2 ^)' /2 1 p _ 25 (2Am) 1/a 1 

V^ 10 / x> o — 2 56 rr~i , , a /m » M 



256 ttVi/ 1+S/T' M>-ff4 ttVv 1+S/T* 

It will be seen later that S is the well-known " Sutb:erland ? s constant" (§11 (F)). 

§ 10. Numerical Calculations for Particular Molecular Models. 

Rigid Elastic Spheres. 

§ 10 (A) In the last section we determined the complete expression for the velocity- 
distribution function for a gas composed of Maxwellian molecules. In the other cases 
there considered we must be content to make numerical approximations, which can, 
of course, be carried to any desired degree of accuracy. We shall consider in most 
detail the case of rigid elastic spherical molecules, for which we shall calculate 
b rs and c rs for = r ■= 3, := s = 3. These are chosen for the fullest treatment partly 
because of their simplicity, and partly as representing the limit between which, and 
the case of Maxwellian molecules, the molecules of actual gases appear to lie. 

In making such numerical approximations the following table of expanded formulae 
for B* (r, s) is useful : — 

Table I. — Expressions for W (r, s). 

B>(0,0)=1 B 1 (l,l)=|^ B 2 (2,2)-^ 5 xy B 3 (3,3) = HsV B 4 (4,4) = ^y 

B° (1, 0) = x 2 + y 2 B 1 (2, 1) = ixy (x 2 + y 2 ) B 2 (3, 2) - %xy (x 2 + y 2 ) 

B 8 (4, 3) = f ixy (x 2 + y 2 ) B 4 (5, 4) = W*Y (x 2 + y 2 ) 

B° (2, 0) - x 4 + ~VW + 2/ 4 B 1 (3, 1) - 2xy ( x * + ^x 2 y 2 + y*) 

B 2 (4, 2) = ^x 2 y 2 (x 4 + ~\*-x 2 y 2 + y*) 

B° (3, 0) - x Q + 7x 4 y 2 + 7xY + y Q B 1 (4, l) - %xy (x Q + \ 7 ~xy + \ 7 ~xY + y 6 ) 

B 2 (5, 2) =^x 2 y 2 (x 6 + ^xy + ^xY + y% 
VOL, cox vi. — A, 2 y 
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It is useful also to recall that B 7c (r, s) = B /c (s, r), and that B A (r, 5) = if 
r <.k or s < k. 

By using these formulae in connection with (143), (144), (147) and (204)-(208) the 
elements b rs9 c rsy have been calculated,^ as above mentioned, and are given in the 
following expressions : — 



(217) V(6„) 



(218) V(c rs ) 



1 


15 


41 


2.7 


2 2 3 2 


15 


269 


5993 


2.7 


2 2 7 2 


2 3 3 2 7 2 


41 


5993 


152537 


2 2 3 2 


2 3 3 2 7 2 


2'3*7 2 


741 


7571 


1517873 


2 3 7.n 

• • • 


2*3.7.11 

a e 


2 5 3 3 7 2 11 


1 


15 


41 


2.7 


2 2 3 2 


15 


877 


6893 


2.7 


2 2 3.7 2 


2 3 3 2 7 2 


41 


6893 


193329 


2 2 3 2 


2 3 3 2 7 2 


2*3*7 2 


741 


3889 


6202777-f 


2 3 7.11 


2 4 3 2 11 


2 5 3 4 7 2 11 


• • • 


see 


»«• 



or, writing out the elements in decimals to six places, 



(219) V(6„) 



1*000,000 
1-071,429 
1138,889 
1*202,922 



1*071,429 

1*372,449 
1*698,696 
2*048,431 



1*138,889 
1*698,696 
2*402,006 
3*259,364 



(220) V (c„) 



1*000,000 
1*071,429 
1*138,889 
1*202,922 



1*071,429 
1*491,497 
1*953,798 

2*455,177 



1*138,889 
1*953,798 
3*044,359 
4*439,790 



741 



2 3 7. 1 1 

7571 
2*3.7.11 

1517873 
2 5 3 3 7 2 1 1 

50375871 
"2 6 3 3 7 2 11 2 



741 
2 3 7.11 

3889 
2 4 3 2 11 

62 02777f - 
2 f, 3*7 2 ll 

225937695 

2 6 3 4 7 2 11 2 



1*202,922 
2*048,431 
3*259,364 
4*916,968 



1*202,922 
2*455,177 
4*439,790 
7*350,929 



* A considerable part of the computations of § 10(A) have been made by Mr. J. Marshall, Scholar 
of Trinity College, Cambridge, who has thus been of much assistance in bringing the results into a 
useful form. 
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As explained in § 8 (E), however, the determinants V (S rs b rs ) and V (S„c rt ) derived 
from the above by successive differencing are much more useful and suitable for 
numerical calculation.* They may conveniently be written as follows, where the 
factor above each column or before' each row is to be multiplied into all the elements 
of that column or row as written (e.g., the right-hand element on the second row of 
V (S„b„) is equal to -59 . {2 3 3.7.11.2.7}- 1 ) :— 



(221) V (SJ)„) 



1 (2.7)" 1 (2 2 3 a 7)- 1 (2 3 3.7.11)- 1 



(2.7)- 1 

(2 2 3 2 7)- 1 

(2 3 3.7.11)- 1 



1 


1 


-1 


1 . . 


1 


45 


103 


-59 . . 


-1 


103 


5657 


6783 . . 


1 


-59 


6783 


149749 . . 



(222) V (<U,) 



(2.7)- 1 (2 a 3 2 7)~ 1 (2 3 3.7.11)" 1 



1 


1 


-1 


1 . . 


1 


2 5 
3 


163 


2,8 7 


-1 


163 


• 11889 


16798f . . 


1 


287 
3 


16798$ 


329573|f • • 



1 

(2.7)- 

(2 2 3 2 7)- 1 

(2 3 3.7.11)" 1 



As we are throughout concerned with ratios of determinants, the above fractional 
expressions for the elements, from which the column -factors or row-factors can for 
many purposes be omitted, are the most suitable for calculation. The following 
values of the elements in decimal notation (to six places) are of interest, however, as 
showing the relative magnitudes of the various terms : — 

(223) V(S rs b rs )= 1-000,000 0-071,429 -0-003,968 

0-071,429 0*229,592 0'029,195 

-0-003,968 0-029,195 0-089,081 

0-000,541 -0-002,280 0'014,565 



0-000,541 
0*002,280 
0-014,565 
0-043,849 



(224) V(3„c„) = 



1-000,000 
0-071,429 
0-003,968 
0-000,541 



0-071,429 
0-348,639 
0-046,202 
0-003,698 



■0-003,968 
0-046,202 
0-187,216 
0-036,072 



0-000,541 
0-003,698 
0-036,072 
0*096,504 



* This process of differencing renders the determinants much more convergent in appearance (cf. (219) 
and (220) with (223) and (224), without really altering in the least their value or the value of any of the 
partial determinants formed by the first n rows and columns. 

2 Y 2 
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By taking these determinants with one, two, three, and four columns we get 
successive approximations to V, V (q/! § 8 (E)), and to the actual coefficients /3 and y, 
as follows : — 



Table II. — -Rigid Elastic Spheres. 



1st approximation 
2nd 
3rd 
4th 



5? 



5J 



>) 



V (b rt ). 



1-000,00 
0-224,49 
0*019,13 
0-000,79 



V (c r8 ). 



1-000,00 
0-343,54 
0-062,15 
0*005,54 



jL/kJq 



V (8 rs b rs ) 

V (8 rs b rs ) ' 



1-000,000 
1-022,727 
1*024,818 
1-025,134 



Zy, 



V (brsCrs) 
V [o rs C rs ) 



1 * 000,000 
1-014,851 
1-015,879 
1*016,065 



The determinants V (8„b rs ), V (S rs c rs ) are obviously much more convergent in form 
than V (b rs ), V (o rs ). Table II. shows that in each case these determinants converge 
rapidly to the value zero, but that the principal minors of the former determinants 
converge also to the same value in nearly constant ratios. These ratios, the 
successive approximations to which are given in the two last columns of Table II. , 
are the quantities X/3 r and 2y r which we require ; they evidently converge rapidly, 
the successive differences being as follows :— 

Table III. — Rigid Elastic Spheres. 





>u kJw** 


Differences. 


/ y* 


Differences. 


1st approximation 


1*000,00 


Aii ( o 


1-000,00 


1485 


2nd „ 


1-022,73 


209 


1*014,85 


103 


OF& 5j 


JL Villjt* \ZiW 


31 


1-015,88 


29 


JEruXJl $« 


1*025,13 




1*016,07 





We may therefore conclude that, within a small fraction per cent., 2/3 r and Ey, 
have the following values for rigid elastic spheres :— 



(225) 



GO 



00 



oo 



00 







1*026, 2y r = 1*016, 2/3 r /2y r = 1*010, 



It should be noticed that even the second approximations to these quantities give 
results which are very nearly accurate, owing to the rapid diminution of the successive 
differences. 
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While we may thus obtain a close approximation to the values of the series 2/3 r 
and 2y r with little difficulty, the approximations to the values of the individual 
coefficients /3 and y converge by no means quickly, as the following table 
shows : — • 

Table IV. — Eigid Elastic Spheres. 





1st approximation. 


2nd approximation. 


3rd approximation. 


4th approximation. 

1*623,0 
-0-943,2 

0-432,8 
-0-087,5 


ft 
ft 
ft 

ft 


1-000,0 


1-340,9 
-0*318,2 


1-520,2 

-0-652,1 

0-156,7 


73 


1-000,0 


1-222,8 
-0-207,9 


,1-309,4 

-0-368,8 

0*075,4 


1-366,3 
-0-526,3 

" 221,8 
-0-045,7 



Evidently the /3's and y's alternate in sign, and successive terms do not seem to 
diminish quickly, at any rate near the beginning. To obtain an accurate estimate 
of the real values of these coefficients it is clearly necessary to carry the approxi- 
mation much further than we have done, but for our purpose this is not required. 



Molecules which are Point Centres of Force varying as r 



— n 



§ 10 (B) The next simplest case, analytically, to that which has just been discussed 
is the case of molecules which are point centres of force varying inversely as the 
n th power of the distance. By comparison of (186) and (205), in conjunction with 
the general expressions for b rs and c rs , it is easy to see that the difference between 

1 

the values of b rs or c rs in the two cases consists of a power series in , the constant 

n — l 

term of which is zero, while the term of highest order is (n— l)~ (w) # Numerically 

the difference is small, as may easily be verified in any particular case ; it appears 

to be of constant sign, b rs and c rs being greatest for molecules which are rigid elastic 

spheres. The behaviour of the determinants V (b rs ), V(c rs ) or V (§ rs b rs ), V (S rs c rs ) is 

similar in the two cases, the convergence being slightly the more rapid in the present 

instance. Since for rigid elastic spheres the second approximation to 2& and 2y r 

proved so satisfactory, we shall be content with a second approximation only, for 

molecules which are point centres of force; this very materially lightens the labour 
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of numerical calculation. The following are the expressions found for the deter- 
minants V (S rs b rs ), V (S rs c rs ) as far as regards the first four elements : — 



(226) V(S„b„) 



1 



1 



7 



14 2(w-l). 



_45 
196 



114 2(w-l). 

4 . -4 



+ 



V (§ rs c rs ) 



1 



49(w-l) 49(n-l) : 



7 



14 



(n-1). 



^~< 



7 



205 



4 



+ 



4 



14 2(w-l)J [588 49(w-l) 49(w-l) 



When w is made infinite these become identical with (221), (222) ; it is interesting 
to notice that the additional terms are the same in the two determinants, though 
whether this is true for other values of r and s is not clear. 

The first approximations to 2/8 r , 2y r are, of course, unity ; the second are found 
to be approximately as follows : — ■ 



(227) 



2/5,. 



45 

44 



1 



16 n-2 



1- 



45 (n-1) 
2 1 



2 






48 n-2 

205 205(n-l) s 

202 



1 



12 



101re-l 

QO 00 

5 the values of 2/9 r and 2y r are 







11 n~~l 

From § 9 (C), (196), we know that when n = 

unity, and this is also true of any approximation to their values made in the present 
manner. From § 10 (A), however, we know that for n = oo the second approxi- 
mations are slightly too small, by 0*003 and 0*001 very nearly. In the following 
table, therefore, which gives the approximate values of 2/3 r and Ey r for various 
values of n lying between 5 and co, the results obtained from (227) have been 
increased by 0*001, 0*002, or 0*003, as seemed most appropriate in each case. 



JL jCv J— # JLj JJJ T • 




/V1JL1U11 diC 1 Ull 


l\j VJCIIIICO UI 


lui uo v cti y mg 


dO 7 




Maxwell's 

case, 

n = 5. 


n = 9. 


n = 15. 


?i = 25. 


ft = OO, 

rigid elastic 
spheres. 


JLmKJn* 




1 


1-007 


1-013 


1-018 


1 • 026 


00 

€ c — _ Zjy r 



1 


1-004 


1-007 


1-011 


1-016 


OO 00 




1 


1-003 


1-006 


1-007 


1-010 



THEORY OF VISCOSITY AND THERMAL CONDUCTION, IN A MONATOMIC GAS. 333 



Rigid Elastic Attracting Spherical Molecules. 

§ 10 (C) Considering next the case of rigid elastic spherical molecules which exert 
attractive forces, as in § 9 (E), we shall neglect the effect of the attractions in 
producing deflections without actual collisions, and, as in § 10 (B), we shall content 
ourselves with a second approximation to 2/3,. and 2y r . The difference between this 
case and that of rigid elastic spheres without attraction is small. The expressions 
for V (§ rs b rs ) V (S r8 c r8 ) 9 as far as regards the first four elements, are as follows : — ■ 



(228) 



V (S n b n ) 



i 1-S/T ] 

T4 1 + S/T 



j 1 — o/ ! 
T? I+S/T. 

45 1 + IS/T ) 
196 1+S/T 



¥"7 /ft \ 

V (c) rs C rs ) 



Ii 1 -S/T \ 

V" 1+S/T J * * ' 



i 1— S/T \ J 20 5 1+ Ml S/T 
^ 1 + S/T J V wg 1+S/T 



When S — 0, i.e., when there is no attraction, these reduce to (221), (223). 
The second approximations to 2/3,, and 2y r are hence found to he as follows :• 

45 ( 1+|S/T) 
= _____^^ 9 approximately, 

205 (] + ^fS/T) 
20^lT^P7^ ' approximately. 



(229) 
(230) 



2/3, 



2 



% 



Since S/T is never negative, it is clear from (229), (230) that the second approxi- 
mations to 2/3 r and 2y r are never less than unity. Their values, without any 
estimated correction for the error of approximation, are given in the following table, 
for various values of S/T. The correction as estimated is appended as a suffix, and 
is to be added to the last digit of the corresponding number. 





X ABLIii V JL." 


— Rigid Elastic Attracting 


y Spherical Molecules 


i. 






Low temperatures. 


Moderate temperatures. 


Hign temperatures. 


s 

T' 


OC. 


5. 

1-015 2 
l-009i 
l-006i 


4. 

1-012 2 
l-007i 
l-005i 


3. 

l*008i 
1*005! 
l-003 


2. 


1. 


0-7. 


0-4. 

1 • 005i 
l-004 
1-OOla 


0-2. 


0-1. 


0. 


CO 
ZjKJ/y 



CO 



GO 00 

o^zz: jdp r f2/y r 




1-038 
1-016 
1-022 


l-004 
1-002 
l-002 


1 
1 
1 


1-001 
l-OOlo 
l-OOOo 


1-Ollg 

l-009i 

l-002i 


1-016 2 
l-013i 
1-003! 


l-023 3 
l-015i 
l-008 2 
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In the case of the actual gases for which S has been determined, it has an extreme 
range from about 50 to 250, while the range of absolute temperature over which 
experiments are usually made is from about 50° C. to 500° C. Thus the limits 5 and -J- 
are rather extreme values of S/T, but from the above table it appears that the variation 
in 2/3 r , 2y r or their quotient hardly exceeds 1 per cent, over this range. The variation 
is especially slow in the neighbourhood of S/T = 1. 

§11. Viscosity and Thermal Conduction. 

We now proceed to apply the expression for the velocity-distribution function 
(§ 8 (B)) to the determination of the coefficients of viscosity and thermal conduction. 
We shall first obtain general formulae for these coefficients, true for any monatomic 
gas, afterwards considering special molecular models in conjunction with the results 
of §§ 9, 10. 

The Coefficient of Viscosity. 
§ 11 (A) The system of pressures at any point of a gas is given by the equation 



(231) P w - pU 2 , P„ - P UV, &c. 

By means of (5) and the velocity-distribution function (142), we find that 



(232) P« = / ,U a = i / ,O ! ' + i /0 3U 8 -O 



2 
? 



00 



2hm I5hm 



— c n G 2y r . 



GO 



(233) V XIJ = - ^JL- c 12 C 2y r . 



Since, by (74), 
we have 



C ll + C 22 + C 33 ~ 0? 



W P „ + P„ +P „ = 3 _£_=„, 

p being the hydrostatic pressure as usually defined. 

By comparing (232), (233) with the equations giving the system of pressures in a 
viscous fluid having a coefficient of viscosity /a, Viz. , with 

(235) P„= 1 ,-i M ( r 2 3ttn dVn ^ 



dx dy dz ) P 



P 3M c n- 



(236) P„ = _„(| + ^) = _ Wu 
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(remembering the meaning of c lu c 12 , &c., as defined in (72)), it appears that the two 
are identical if we write 



00 



(237) M = Tt)7 — C 2y r . 

Hence, according to the kinetic theory, a gas behaves like an ordinary viscous 
fluid having a coefficient of viscosity defined in terms of the molecular data by (237). 
By (170), (173) we have 

(238) C = 25 
whence, also, 



TTVK, 







5 1 °° 

(239) p. = — — • 2y r . 

JjIT JIKq 



CO 



As we have seen in § 8 (B), k and Sy r are functions of the temperature (or h) only, 



and v does not appear at all in the formula for /a. Hence, within the limits of 
applicability of our theory (cf § 2), the coefficient of viscosity of a gas is independent 
of its density, varying only with the temperature. The law of this variation depends 
on the law of inter-action between two molecules at encounter, this being involved 
through cp k (ry). As this function has remained unspecified, the expression in f237) 
is perfectly general and valid for any nearly perfect monatomic gas. 

The Equation of Energy for a Simple Monatomic Gas. 

§ 11 (B) In the discussion of the equation of transfer in § 3, we consistently neglected 
such second order quantities as products of differentials, or differentials of small 
quantities like UV, C 2 — 3U 2 , and so on. In this way we have obtained an expression 
for the velocity-distribution function which is correct to the first order. By means 
of this function we can now determine the values of UV, C 2 — 3U 2 , UC 2 , and similar 
expressions which are of the first order of small quantities, and by substitution in the 
equation of transfer obtain this in a form accurate to the secpnd order. This we shall 
do for the special case Q = (u) 2 + (v) 2 + (w) 2 , in order to get a second approximation to 
the equation of energy. 

From the velocity-distribution function, using the formula (237) for the coefficient 
of viscosity, we have 

(240) C 2 - 3 (hm)- 1 = 3E,T/m. 

du 



(241) 9U 2 -C 2 =-2{rf P )(2 



dx 



(242) OV — (M (fa + dy 

VOL. COX VI. A. 2 Z 



3v . 


dw a 


dy 


dz, 


du \ 
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In finding U C 2 it is convenient to write 



(243) 

and to recall (cf. (170)) that 

(244) 

Thus we have 

(245) 



00 



00 










-L>0 — " 4M)* 



1 3T "r + 4 



UC 2 = -f Bo {2hm)- 2 ± -^ 2 '-J-* /3 r _ 



r 3a; r 

-iB (2M- 2 ^-m 

1 dx o 



¥ 



[x R 3T 
p m dx ' 



where we have eliminated /3_i by means of (148). 
Again, if Q = (uf + (vf -f (wf , we have 



Q 
Q 



«o 2 + V + ^o 2 + 2 ( Wo U + v V + w W) + C 2 , 



!i .J, ...8, ni 



(w)Q 

3< 



3(w) 



w +V+w +C'. 

w (m 2 + V + w 2 ) + u U 2 + 2 (> U 2 + v UV + w UW) + UC 2 

(w), 



3Q 



d (w), 



6vAe 



Hence, putting w = v = w = except in differential coefficients, we have 



_3 
dt 



(»Q) 



v 



dv 3 

T" OP 



L\ 



dt \2hmJ 2hm dt l dt \2hmJ 



o R^T fdu f 3t> f 3m 






m \3x 3^/ 






82; 



+ 3-' 



Ki/ 3T 



l x (Au)Q) 






a 3^o 



.2 C 2 +2U^+2DV 



3t? 



m 3x * 



3 m- 



* 3 






3x 



12 3^0 



= ,2 fC» + f3U 2 -C 2 ^ + 2UV^ + 2uW 



^ + 2UW~^° +2^(,UC 2 ) 

ox dx dx 



du , nT — dv n . ^T-rrx-r 3m, 



3 



3x 
R*/T / du a . dv 



dx 



dx 



. . +2^-(,UC 2 ) 

dx J 3x 



o i ^^o , 3m 



m \ 3a? 3?/ 3z / 



P 



Ay { ° U A 4 ' ^ 
1 3aj/ ~"" 3 



3^ 
3x 



2 



'dv , 3w Vl , ^ 3 



^+^p +2f (,U0 2 ). 
9z 3w/ J 3ic 



3 (-u), 



0. 



Also, since no energy is gained or lost in molecular encounter 



s. 



AQ = 0. 
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The equation of transfer becomes, consequently, after a little reduction, 

(oAa\ 3 T _ v 9 //. 3T\ 2 rp/duo dv n dw \ 

(246) P - - 2 -(/ M ^) -spT{-+~ + — 



2 m^ f 9T /3w \ 2 g. / v 3w o y y /S-Wo 3wc ' 

3 R7 r Z W "H ¥j "^Vfo + "ty 



which is the equation of energy. 

The Thermal Conductivity of a Gas. 

§ 11 (C) In the equation of energy which we have just obtained, the second term 
on the right-hand side represents the change of heat per unit volume due to the 
variation in density at the point considered, while the third term may be proved 
equal to the heat produced by internal friction. The first term, by comparison with 
Fourier's equation of conduction of heat (S- being the thermal conductivity and G v the 
specific heat at constant volume), i.e., with 

P 3T _ y d^ /q 9T v 

v dt dx\ dx, 

is seen to represent the change of heat by conduction, and to indicate that the 
coefficient of thermal conductivity of a gas is given by 

(247) &=/a*C.. 

The value of / in this well-known formula is, for a general monatomic gas, given by 
(243), i.e., 

oo / oo 

(248) /=|2A./2 yr . 



0/0 



In general f is a function of the temperature only. 

Formulce for /m and $■ for Particular Molecular Models. 

§ 11 (D) By substitution of the values of C and 2/3 r /2y r given in §§ 9, 10, for the 
particular molecular models there discussed, we obtain the following special cases of 
(237) and (248) :— 

(249) Rigid elastic spheres, 

M = roi6 ^^J-^T f= f. l'OlO = 2*525, 

(250) Attracting spheres, 

M = (1+ea) 64&w T ) rrs/T y = *(i+*.)> 

2 z 2 
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(251) Centres of force r~ n , 

M = (1 + 75m (R T ^ /= |(l + 4). 

\ n— 1, 
In these formulae e a and e c denote the values of 2y r in the cases respectively when 



the molecules are attracting spheres and centres of force, and S a and S c similarly denote 
2/3 r /2y r ; their values are given in Tables V. and VI., and in no case differ from 
unity by more than two per cent. 

The mode of variation of jul with the temperature affords a guide to the law of inter- 
action between the molecules of actual gases. By comparison with experimental 
determinations of jul at various temperatures it is thus found that of the above models 
the one which most closely represents the behaviour of actual molecules in this respect, 
at ordinary temperatures, is the second, i.e., a rigid elastic attracting sphere. # 

Comparison of the present formula for jul and § with those of my previous paper. 
§ 11 (E) The general formulae (237) and (248) for viscosity and thermal conductivity 

00 CO 00 

agree with those of my former paper, t except that the factors 2j8 r and 2/3 r /2y r were 



there omitted. This was in consequence of the assumption on which the analysis of 
that investigation was based, that F(U, V, W) is sufficiently represented by the 
terms of the first three degrees in U, V, W. We have seen in § 9 (C) that this is true 
for a gas composed of Maxwellian molecules, but not otherwise. It seems of interest 
to consider why the % neglect of all the coefficients /3 r , y r after r = led to results of 
such accuracy ; for the errors arising from the assumption are represented in the 
special cases (249)-(25l) of (237) and (248) by the factors 1*016, l'OlO, l + € fl , 1 + 3 M 
1 -he c , 1 + <J C , so that the necessary corrections to my previous formulae do not exceed one 
or two per cent. Enskog-, on the other hand, after deducing formulae similar to (237), 
(248), but without evaluating the coefficients /3 r9 y r , made a first approximation by 
neglect of all these coefficients after r — 0, and arrived at the result f= 5 for rigid 
elastic spheres. J This was due to the fact that such a use of (237), (248), as they 
stand, involves not only the neglect of all the coefficients after r = 0, but also requires an 
assumption as to the values of /3 , y themselves, as, for instance, that they are approxi- 
mately the same as for Maxwellian molecules ; a comparison of (196) with Table IV. 
(p. 331) will show that this is far from being the case. 

It may readily be seen, however, that the method of my former paper required no 

* At very low temperatures, however, the n th power centre of force is the molecular model which gives 
by far the best representation of the relation between //, and T, in the case of helium ; cf. Kamerlingh 
Onnes and Sophtjs Weber, 'Comm. Phys. Lab. Leyden/ 134b, p. 18, or Jeans' < Dynamical Theory 
of Gases/ 2nd ed., §§ 405, 407. 

t Chapman, 'Phil. Trans./ A, vol. 211, p. 433, et seq. 

t Enskog; ' Phys. Zeit./ XII, p. 58, 1911. 
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hypothesis as to the values of /3 , y , which were determined from AUC 2 and AU 2 as in 
this paper ; in effect, the later coefficients were neglected, while fi b and y were obtained 
from the equations (138), (139) corresponding to s = 0. Viewed otherwise, my 

previous formulae were equivalent to (237), (248) divided respectively by 2/c r y r //c ( 








a> oo 



and 2/c r /3 r /2/c r y r . Thus the neglected factors in jm and f were 



OO CO oo 

lE,K y r 2/3 r ^£;K r y r 

£ and £ — £ , 

CO OO OO ' 

22/c r y r 2y r 2/c r p r . 

o o ' o 

which evidently reduce to unity if we neglect all the /3' and y\ after r = 0, without 
any assumption as to the values of /3 and y . 

One of the main results of the former paper was that f \= -§ for all monatomic gases, 
and not only for those composed of Maxwellian molecules. This is now seen to require 
modification, but the values here found for f in the special cases which have been 
considered in §§9, 10 show that the correction needed to make the equation accurate 
is very small ; it appears probable that for all likely molecular models f is very 
slightly greater than 2*5, and that it is nearly but not quite constant with change 
of temperature (except when the molecules are elastic spheres or centres of force 
proportional to r~ n ). 

Comparison of the Formulae for m and S- with the Results of other Theories. 

§ 1 1 (F) The only kinetic theory of viscosity and thermal conductivity which could 
hitherto lay claim to numerical accuracy (within the limits imposed by the initial 
postulates) is Maxwell's theory* of a gas composed of molecules of the kind dealt 
with in § 9 (C). The results of his theory are special cases of the general formulae of 
this paper. 

The theory of a gas composed of molecules which are point centres of force varying 
inversely as the n th power of the distance had not been discussed in detail, previous 
to my own former paper. B,AYLEiGH,f however, from considerations of dimensions 
alone, had deduced the law of variation of viscosity with temperature, and the same 
argument would also show that for such a gas / is an absolute constant (for any given 
value of n). Nothing was known as to the value of this constant, or of the numerical 
coefficient in the expression for /*, and it is a surprising result, which could hardly 
have been guessed a priori \ that as n ranges from 5 to oo the value of/ should vary 
only from 2*500 to 2*525 approximately. 

The theory for molecules which are rigid elastic spheres exerting attractive forces 
was equally undeveloped. Suthekland J had taken an important step, however, in 

* Maxwell, 'Collected Papers/ vol. II, p. 23. 
t Eayleigh, 'Boy. Soc. Proc.,' vol. 6, p. 68, 1900. 
X Sutherland, < Phil. Mag.,' (5), 31, 1893. 
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deducing the correction to the law connecting fx and T (i.e., jm oc T 1/2 ) for molecules 
which are rigid elastic spheres without attraction ; he showed, without attempting 
numerical accuracy, but by a method which is correct in its main outlines, that the 
attractive forces necessitate an additional factor (l + S/T)" 1 , as in (249). The law 
jii oc T 1/s (l +S/T)" 1 is more successful than any other in representing the observed 
relation between jm and T over a considerable range of temperature, and S is deservedly 
known as Sutherland's constant. 

The theory of a gas composed of molecules which are rigid elastic spheres, which 
was taken by Sutherland as the basis of his modified formulae, has been developed 
along lines different from those of this paper by Clausius, Maxwell, Boltzmann, 
Meyer, Stefan, Jeans, and others. Their method was less analytical than the 
present one, and while it gave correctly the general relationships between /u, 9-, p, p, 
and T, its results do not possess numerical precision. Jeans* notably improved 
certain of the formula due to earlier authors by taking into account the tendency 
of a molecule to persist, after a collision, in the general direction of its original course. 
For this reason his expression for the viscosity, viz., 

(252) M = _ — — ( — I , (Jeans) 

4x ir h cr\m j 

approaches more nearly to the correct expression (249) than does the formula of any 
other author t. A comparison of (249) with (252) indicates that the latter is still too 
small by 12 per cent. ; the error of the original formula, without Jeans' correction, 
was 30 per cent. 

The numerical inaccuracy of the earlier prevailing theory of conductivity, which 
was due to Meyer J, was very great. Its result was generally given as 

& =//aC v where /= 1*6027, 

but Prof. L. V. King, of McGill University, has pointed out to me by letter that 
Meyer's argument really leads to the result 

jf= 1*4161, 

a numerical mistake having crept into his work which had not previously been 
detected. The correct value of/' for rigid elastic spheres is given in (249), i.e., 

/= 2'525. 

This large error in Meyer's theory indicates the difficulty of arriving at numeri- 
cally accurate formulae by the older " mean free path " method, and diminishes 

* Cf. Jeans' ' Dynamical Theory of Gases.' 

t Apart from that in my former paper, which was 1 * 6 per cent, too small. 

J Meyer's ' Kinetic Theory of Gases/ 2nd English edition, chap. IX. 
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confidence in other parts of that theory where detection of error is less easy. Until 
recently Meyer's value of f received support from experimental data for diatomic 
gases, to which it does not really apply ; only lately have data for monatomic gases 
been obtained, which, as we shall see, give values of/* nearly equal to f. 

§ 12. Comparison of the Theory with Experimental Data. 

The Variation of Viscosity with Pressure, 

§ 12 (A) The main objects of a comparison of a mathematical theory with experi- 
mental data are either to test whether the postulates underlying the theory are valid, or 
whether the theory is itself mathematically correct. The present theory being exact, 
within certain defined limits, our purpose in this chapter is to consider how far the 
hypotheses underlying the analysis are well founded. The general validity of the 
foundations of the kinetic theory is attested in many ways, one of the most striking 
being the independence of viscosity and pressure in a gas. This law, when first 
discovered by Maxwell, seemed so improbable that it gave a great stimulus to 
experimental research on gases, and the constancy of /a, when T is kept constant, has 
been verified over a range of pressure extending from a few millimetres of mercury up 
to more than one atmosphere. Warburg and von Babo have found that, in the case 
of carbon dioxide, the law begins to fail when the pressure becomes so great as 30 to 
120 atmospheres, /x rising appreciably. In very rarefied gases, on the other hand, the 
viscosity falls below the value appropriate to the existing temperature. This must be 
referred to the failure of the postulates of our theory to represent the facts in these 
extreme cases, the molecules becoming too few for our statistical method to apply, on 
the one hand, while on the other our assumption that the molecular paths are 
rectilinear for the major part of the time, and our neglect of multiple encounters, 
become illegitimate. 

The Variation of Viscosity with Temperature* 

§ 12 (B) Over a wide range of pressure and temperature, undoubtedly, the general 
postulates of our theory are true for actual gases. We cannot discover directly, 
however, the nature of the molecules or their mode of collision, and it is important, 
therefore, to examine which molecular model yields formulae most in accordance with 
experimental data. For this purpose we naturally choose those properties which are 
most affected by the nature of the molecule ; the chief of these is the variation of 
viscosity with temperature. Maxwell abandoned his theory of a gas composed of 
rigid elastic spherical molecules because it led to the relation jul oc T 1/a , while his experi- 
ments gave the result /x oc T. This caused him to develop the theory of a Maxwellian 
gas (§ 9 (C)), for which jul cc T, but later experimenters have failed to confirm this law, 
and we must conclude that the molecules of actual gases behave during encounters 
neither like elastic spheres nor like Maxwellian molecules. The observed relation 

* The reader may be referred with advantage to the discussion of this point by Jeans in the second 
edition of his ' Dynamical Theory of Gases/ §§ 399-407. 
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between fx and T agrees much more closely with Sutherland's law jul oc T^l + S/T)" 1 
than with any other ; for example, Schmitt^ has found that the law is valid for 
hydrogen and helium from —60° C. to 185° C, and' Barus has shown that it holds 
good for air over a wide range of temperature. The law has not been tested, for the 
former gases, above 185° C. 

This seems to indicate that for the kinetic theory of gases at ordinary temperatures 
the best molecular model is an attracting sphere, and it is interesting to notice that 
this model is the one used by van der Waals with such success in deducing his 
famous law. Further confirmation is supplied by the excellent agreement between 
the values of the molecular diameters deduced on this hypothesis from the constant 
b of van der Waals' law and from the viscosity by means of my formula (250) — ■ 
of- § 12 (F). 

At low temperatures Schmitt^, BestelmeyeivI", VogelJ and others have shown 
that the observed values of /x are greater than those predicted by Sutherland's 
law. This may be compared with the rise in the value of jul when the pressure is 
greatly increased, both effects probably having a like cause ; in these states, when the 
mean free path of the molecule is much reduced, the molecular paths may cease to be 
approximately rectilinear between collisions, and multiple encounters will grow in 
importance. Since our theory rules out these contingencies, its results cease to be 
applicable, and a modification of the theory and its postulates is necessary if a proper 
account of these phenomena is to be given. In regard to this, one point which should be 
noticed is that in § 9 (E) a term f 2k (y) in <j> 2k (ry) was neglected (of. (211)) which, if 
retained, would cause the law connecting jul and T to take the form 

T J/a 

/xoc 



1+(S/T)+/(T) 

where f (T) can be expanded in the form AT~ 2 + BT~ 3 + .... This term is due to the 
effect of the attractive forces in producing deflections without the occurrence of 
collisions, and is probably always small ; but it may readily be seen that it is always 
positive, and that this correction would lead to &■ diminution in the theoretical value 
of jul at low temperatures. Clearly, therefore, the observed discrepancies cannot be 
attributed to our neglect of this small quantity. § 

* Schmitt, 'Ann. d. Phys.,' 30, p. 399, 1909. 

t Bestelmeyer, 'Munich dissertation/ 1903. 

J Vogel, 'Berlin dissertation,' 1914, where full references, and an interesting discussion of low 
temperature work on viscosity, are given. 

§ Vogel, in his dissertation, suggests as possible causes of the failure of the theory to represent the 
observed variation of jjl with T at low temperatures (i.) a failure of the ordinary mechanics, such as is 
contemplated in Planck's theory of quanta ; (ii.) that the attracting sphere model no longer represents 
the molecule ; (iii.) that 1 + S/T should be replaced (according to my suggestion in ' Phil. Trans./ A, 
vol. 211, p. 474, 1912) by 1 + (S/T) ± (C'/T) 2 . By the latter means a better accordance with observation 
is obtained, but the new term has the minus sign, and is therefore illegitimate. 
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The Thermal Conductivity of Monatomic Gases. 

§ 12 (C) It is convenient to discuss the thermal conductivity of gases in terms of the 
constant f in the formula 3- —//ulC v , as this eliminates the necessity for a separate 
discussion of the dependence of 9- on pressure and temperature ; this is parallel with 
that of ya, and f is nearly or quite independent of pressure and temperature in normal 
conditions. As we have seen in § 11 (F), the value of f has been a matter of some 
uncertainty ; so long as its value for rigid elastic spheres was supposed to be 1*6027, 
while for Maxwellian molecules it was known to be f , it seemed to offer a means of 
testing the suitability of different molecular models. On the ground of the 
discrepancy between the theoretical and observed relation between jjl and T, 
Maxwellian molecules were known to be unsatisfactory representations of actual 
molecules. Until about 1900 no reliable determinations off had been made for mona- 
tomic gases, and those found for diatomic gases agreed fairly well with Meyer's value 
of f (i.e., 1*6027 or, more accurately, 1*416); at the time this was regarded as a 
confirmation of the rigid elastic spherical model of the molecule, and as indicating 
that the internal molecular energy, which is not taken into account in these theories of 
a monatomic gas, is transmitted at the same rate as the translational energy. When, 
in 1902, Schwarze obtained the values of/* for argon and helium, and found them 
nearly equal to f , the conclusion to be drawn was not obvious. It certainly 
contradicted Meyer's theory, but left the question open as to whether the analysis, or 
the assumption of the rigid elastic spherical model, was at fault ; also iff = f indicated 
that the molecules are Maxwellian, the failure of the corresponding law connecting 
juL and T remained unexplained. It should be remembered, moreover, that the law 
/xocT 1/s for rigid elastic spherical molecules is equally contradictory to experiment. 

These difficulties were removed by the theorem of my former paper, according to 
which f is an invariable constant f for all monatomic molecules. This is now seen to 
be incorrect as a general theorem, but the deviations found for the various particular 
molecular models discussed leaves little room for doubt that /is very nearly equal to 
f in the case of all likely models. The fact simply is, therefore, that f is very 
unsuitable as a means of discrimination between different models, and Schwarze's 
observations indicate some mathematical fallacy in Meyer's theory, without supporting 
any particular molecular model. The observed values of f are hardly known with 
sufficient accuracy to enable any conclusion to be drawn from a slight divergence 
from the value |, within the limits prescribed in (249) to (251). They are important, 
however, as confirming the general validity of the kinetic theory, apart from any 
hypothesis as to the nature of the molecules. 

The following table contains all the available data concerning the value off for 
monatomic gases. Only very recently has the conductivity of neon been deter- 
mined, owing to the scarcity of the gas ; for krypton and xenon its value is still 
unknown. 

VOL. COXVI. A. 3 A 
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Table VII. — Values of /for Monatomic Gases. 



vJclS. 



Absolute temperature. 



Helium 

Argon . 

Neon . 



°C. 

273 

81 
21 

273 

91 

283 



{ 



t 


/• 


Authority. 


2-51 


SCHWARZE* 


2 


•40 


EUCKENJ 


2 


'23 


EUCKEN 


2 


•02 


EUCKEN 


2 


50 


SCHWARZE 


2 


49 


EUCKEN 


2- 


57 


EUCKEN 


2-50 


DornJ 



These results for argon and neon and, to a less extent, for helium at normal tempe- 
ratures agree very well with the theoretical value of jT, especially since the combined 
experimental errors in their determination may easily exceed one per cent, at ordinary 
temperatures, and much more at low temperatures. 

The diminution in the value of f for helium at low temperatures, if confirmed by 
further experiment, is very interesting and important. Helium is peculiar at low 
temperatures also in the striking divergence of its viscosity from Sutherland's law. 
Eucken suggests as the explanation of the former phenomenon a partial failure of 
interchange of molecular energy at collision, but (cf. Table VI. of his paper) down to 
81° C, at any rate, the value of G v for helium remains constant and appropriate to a gas 
which possesses only translational energy. A failure in interchange of translational 
energy would contradict the ordinary dynamical laws, and it is certainly desirable to 
seek some other explanation, if this be at all possible. 

The alternatives are not numerous, and will be examined in turn. We may rule 
out a numerical error in the theory, of more than one per cent., as being quite 
improbable ; but though all the molecular models discussed in this paper lead to 
values of f equal to or slightly greater than 2' 5, it is conceivable that for some 
peculiar model f may have rather different values and a wider temperature range. I 
think this is unlikely, and that it is probably possible to prove that f always 
exceeds 2'5, but this is only a speculation ; helium agrees so well at high tempera- 
tures, however, with Sutherland's law connecting /a and T, that its molecules can 
hardly be supposed so different in behaviour from, rigid elastic attracting spheres as 
to make f theoretically equal to 2*0 at low temperatures. 

Again, molecular aggregation might seem to afford an explanation, since if part of 
the gas were polyatomic through clustering of the molecules, the value of f would 

* Schwarze, 'Halle dissertation, Ann. d. Phys./ (4), 11, p. 303, 1903. 

f Eucken, 'Phys. Zeit.,' 14, p. 324, 1913, Tables 3, 6. Eucken (footnote 4 to p. 328) states that 
Sciiwarze's value of /for helium is too large owing to a miscalculation in determining C v . 

J This result was kindly communicated to me by Prof. Dorn, of Halle, as an extract from ' Mitt. d. 
Naturf. Ges. z. Halle,' 4, 1914. 



THEOEY OF VISCOSITY AND THERMAL CONDUCTION, IN A MONATOMIC GAS. 345 

probably lie between 2*5 and the lower values characteristic of polyatomic gases. In 
this case, however, the value of G vi the specific heat at constant volume, should rise 
to correspond with the internal energy of such molecules ; as the experiments indicate 
a constant value of G v , the suggestion must be abandoned. 

The only possible remaining hypothesis seems to be to attribute the fall in /to the 
neglect of multiple collisions between molecules, including also the effect of the 
attractive forces (in Sutherland's case) in producing deflections without collisions ; 
at low temperatures the molecules may be too close together for these postulates of 
our theory to continue valid. If we determine 9- for helium from the formula 2' 5 jul C^, 
using the value of jul calculated from Sutherland's formula (which is less than the 
observed value at low temperatures, as we have seen), the result is less than that 
observed at low temperatures. Hence both 9- and jul diminish with temperature less than 
is predicted by Sutherland's law, the divergence being greater for jm than for S-, so that 
f also diminishes. We cannot enter here into a test, by calculation, of this suggested 
hypothesis, but some confirmation might be sought experimentally by examining 
whether f is less than 2*5 for helium at normal temperatures but under considerably 
increased pressure. The latter would bring the molecules closer together in the same 
way as would a diminution of temperature, and this is all that our suggestion requires. 
It is known that over a large range of pressure jm and & are constant, but that at 
high pressures jul increases ; the independence of & on pressure has usually been tested 
by diminishing the normal pressures, # and experiments under increased pressure 
might throw valuable light on the present phenomenon. Gases other than helium 
may be expected to behave similarly, though perhaps only with lower temperatures 
or higher pressures. 

§ 12 (D) The case of mercury vapour may also be mentioned, as it was the first 
monatomic gas for which/ was determined. Kocnt determined jul for mercury vapour 
at 203° C, 273° C, and 380° C, while SghleiermagherJ determined 3- at 203° 0. 
These data, together with the theoretically calculated value of G v9 led tof— 3*15. 
Meyer and others have raised objections to the determinations of /x (a) because the 
three values show an improbable amount of variation with temperature, and 
(b) because of the vitiating effect of condensed mercury on the walls of the capillary 
tube used in the experiments. Vogel§ has made a re-calculation of jul for mercury from 
an interesting formula which he gives, and finds that at 573° C. absolute || jul should 
equal 593'10~ 7 ; this, combined with Sohleierm acker's result, reduces/ to 2*80. But 
it is desirable that more accurate experiments should be made in order that a 
thoroughly reliable value off may be obtained. 

* Eucken, 'Phys. Zeit.,' 12, p. 1103, 1911, Table 2. 

t Koch, ' Wied. Ann./ 19, p. 857 (1883). 

J Schleiermacher, 'Wied. Ann.,' 36, p. 346 (1889). 

§ Vogel, 'Berlin dissertation/ p. 57, 1914. 

|| So given by Vogel ; it may be a misprint for 473° G. 

3 A 2 
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The Thermal Conductivity of Polyatomic Gases at Loiv Temperatures. 

' § 12 (E) The formula 3-=/^^ is true for polyatomic as well as for monatomic gases, 
f being independent of pressure and temperature over a considerable range. Under 
normal conditions, however, its value is 2*0 or less, being greatest for diatomic gases 
and diminishing down to about 1*5 for complex molecules. EuOKEisr # has made the 
interesting and important discovery, however, that diatomic gases show an increase 
in f at low temperatures, the conductivity varying with temperature in the sense 
opposite to that observed in the case of helium. This is apparent, to a slight extent, 
in nitrogen, but is most striking in the case of hydrogen. It is found that, simul- 
taneously with the rise in f the specific heat G v progressively falls in value until at 
21° C. absolute its amount is that appropriate to a monatomic gas of the same molecular 
weight. At these low temperatures the rotatory motion of the molecules seems to 
fail, for some reason as yet undiscovered, so that the gas behaves in certain respects 
as if its molecules were of the spherically symmetrical type discussed in this paper. 
It is highly interesting and significant that this approach to monatomicity is 
accompanied by an upward tendency of f towards the value (2'5 approximately) 
which is appropriate to monatomic gases. The same phenomenon may be expected in 
the case of the other diatomic gases, at lower temperatures corresponding to their 
lower boiling points. In the following table # the results for hydrogen alone are 
given ; the number n in the third column represents the number of " degrees of 
freedom " of the molecule, as calculated from the observed values of CL. 



Values for / for Hydrogen. 



Absolute 
temperature. 


/. 


n. 


273 

81 
21 


1-96 
2-09 

2-37 


4*80 
4-41 
3-16 

2-98 



The Diameter of the Molecule. 

§ 12 (F) In my former paper tables were given showing the values of the molecular 
diameters for several gases, calculated on the hypothesis that the molecules are rigid 
spheres, with or without attractive force. These require a small correction to be 
strictly accurate, on account of tkfe factor (2y,.) 1/a there omitted from the formula for 

* Cf. Etjcken, <Phys. Zeit./ 14, p. 329, 1913, Table 6. 
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o- in terms of jul. In this paper T shall give only a few values of a-, calculated from 
recent data and from the exact formula 



(i+O 



5m 



/R T v. 



64tt 1/ > (1 + S/T)\m 

0*491 (1 + QpC 
2 1 W(1+S/T) 

appropriate to attracting spherical molecules. The values of S are determined from 
the variation of viscosity with temperature, G is the mean molecular velocity, while 
e a is found from Table VI. (in all the cases here considered it is quite negligible). 

In the following table values of the diameters calculated from the constant b of 
van dee, Waals' law are also given for comparison. The agreement between the two 
sets of values is in most cases remarkable, and the table as a whole is a testimony to 
the close numerical accuracy now attained by the kinetic theory ; where there is 
disagreement in the table there is in most cases uncertainty as to the data. 

While exact agreement may be expected only for monatomic gases, the values for 
diatomic gases show that our theory gives a mean diameter, in other cases, which agrees 
with that found for b ; the internal energy which prevents the application of our 
formulae to the conductivity of polyatomic gases hardly affects viscosity. 

Table VIII. — Molecular Diameters Calculated from Viscosity and 

Van der Waals' Law. 



Gas. 



Argon . 
Krypton 
Xenon . 
Helium . 
Oxygen 
Hydrogen 

Nitrogen 

Air . . 

Carbon dioxide 



/vio 7 . 


S. 


2107 


162 


2334 


182 


2107 


252 


1885 


75 


1923 


130 


854 


76 


1672 


112 


1721 


111 


1388 


249 



l + e a . 


(Table Y.) 


1-002 


1-001 


1-000 


1-006 


1-005 


1-006 


1-003 


1 • 004 


rooo 



Van der Waals' 

b. 















001347 
001774 
002304 
000432 

00096 
00255 
00209 
00228 



2<Ti 


•10 8 


(viscosity). 


2 


•84 


3 


■12 


3 


•47 


1 


•89 


2' 


93 


2 


•36 


3 


10 


3 


•08 


CO CO 


•27 


•20 



2cr 2 -10 8 
(Van der Waals.) 



2 
3 
3 
1 

2 

3 
3 
3 



85 

14 

42 

96 

89 

52] 

54 

08 

30 

40 

20 



Eeferenoes. 

Viscosity jn at 273° C absolute. — These values are taken from the table on p. 476 of 
my first memoir, where full references may be found. They agree very well with the 
list given by Eucken (' Phys. Zeitschr./ 14, Table 3, 1913), in which Vogel's 
determinations are included with other recent values in taking means. 
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Sutherland's constant S. — These values are those given in the same table of my 
former memoir (what is called S in this paper was there, and is usually, denoted by C), 
where references to sources may be found. The value for krypton was there given 
incorrectly ; I am indebted to Dr. G. Rtjdorf for the correction. 

van der Waals' b. — The values for helium (Kamerlingh Onnes), hydrogen, 
nitrogen, air (Eose-Innes), and carbon-dioxide (van der Waals') are taken from 
Jeans' 'Dynamical Theory of Gases,' 2nd edition, § 194. Those for argon, krypton, 
and xenon (Kamsay and Travers), from Kudorf, ' Phil. Mag./ June, 1909, 
p. 795, are not direct experimental values, however, but are calculated from critical 
data. 

Diameter 2^ from viscosity. — These are practically twice the values for the radii 
given on pp. 476, 481 of my former memoir, where, however, errors of calculation 
(here corrected) were made in the case of argon and krypton (as Dr. G. Rtjdorf kindly 
indicated to me). 

Diameter 2<r 2 from van der Waals' b. — The value for oxygen, and the second 
values for nitrogen and carbon-dioxide (as well as 2^ for the latter) are from Table 7 
of Eucken's paper ; he does not give his authorities, but his values are probably the 
most recent and reliable. The value for hydrogen, he says, is doubtful. 



